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FLOOD PROTECTION OF CANALS BY LATERAL SPILLWAYS 


Harald Tults,! A.M. ASCE 
(Proc. Paper 1077) 


SYNOPSIS 


To avoid topping of the banks of canals during flood, lateral spillways are 
‘requently used. However, a lateral spillway by itself is not capable of re- 
jucing the height of the excessive energy head. The energy reduction has to 
9e accomplished upstream from the spillway. 

It has been proven that the water level in front of a lateral spillway crest 
san be computed by the Bernouilli and the continuity equations. No trials are 
required in computing the action of a lateral spillway of uniform section, 
when the friction effect is negligibly small, or the bottom slope is made equal 
0 the average friction slope. 

A brief survey is given about different water surface profiles upstream, in 
ront, and downstream from the spillway, depending on the flow stage and on 
he location of the spillway in the canal. 

The effect of non-uniform velocity distribution is taken into account when 
-omputing the surface profiles in front of the spillway crest. 

To demonstrate the effectiveness of lateral spillways, several spillway 
ayouts with different sections and with additional orifices are analyzed. 
[heir performances are illustrated by depth-flow curves. The conditions 
maintaining the subcritical flow stage in lateral spillway section are dis- 
‘ussed. Further, the efficiency of the combined action of the intake and 
ateral spillway in raising the headwater during floods is investigated. 


Notation 


The letter symbols introduced in this paper are defined where they first 
ppear, in the text, or by illustration, and are assembled alphabetically in the 


Appendix for convenience of reference. 
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A. Hydraulics of Lateral Spillway 


1. Introduction 


Using an open canal for diversion from a reservoir, the canal depth re- 
quired for normal flow has to be increased so that it may convey the maxi- 
mum inflow during the floodwater without overflowing its banks. 

It is economical to keep this freebord of a diversion canal as low as pos- 
sible, particularly, if the canal is of substantial length. This is possible by 
using a lateral spillway which provides the outlet for surplus water. 

Even, when the diversion inflow is regulated by a float or remote con- 
trolled gate, the installation of a lateral spillway is necessary. Such a con- 
trol is not absolutely reliable because of the possible mechanical or power 
failure. 

Further, a lateral spillway is necessary at the end of a canal to discharg 
the full flow at the abrupt closure of the turbines and along a long canal to 
spill the excess inflow from the hillside. 

The first publications of the laboratory tests about the performance of 
lateral spillways by H. Engels2 and G. S. Coleman and D. Smith3 were rathe 
confusing for hydraulic engineers. The tests by Engels indicated a rising 
and those by Coleman and Smith a dropping water surface along the spillway 
crest. This contradiction was caused by the different stages of the flow use 
subcritical in Engels’ and supercritical in Coleman-Smith’s tests. | 

Ph. Forchheimer4 considered the energy level to be parallel to the spill- 
way crest and to the canal bottom and assumed the water surface profile 
along the spillwa ay crest to be a straight line. 

G. De Marchi?” proved theoretically that the energy head along the lateral 
spillway is essentially constant and the water surface profile curved: rising 
in subcritical and dropping in supercritical flow stages. Further, De Marchi 
introduced a graphical method for designing lateral spillways. 

Br. Gentilini6 verified De Marchi’s theory experimentally. 

Jos. Frank’s” publication concerned with the reciprocal action of the in- 
take gate and the lateral s alter. 

O. Herz8 and H. Favre9,10 established general methods for stepwise com 
putations of the water surface in the lateral spillways. 


2. “Mitteilungen aus dem Dresdner Flussbaulaboratorium,” a/Forschungsal 
beiten des Vereins Deutscher Ing., Vol. 200/201, 1917. b/Zeitschrift des 
Vereins Deutscher Ingenieure, 1918, pp. 362, 387, and 412, 1920, p. 101. 

3. “Discharging-Capacity of Side Weirs.” Institution of Civil Engineers, Vc 
No. 6, 1923, London. 

4. “Grundriss der Hydraulik,” 1926, p. 406, Leipzig-Berlin. 

5. “Saggio di teoria del funzionamente degli stramazzi laterali (o stlaratam’ 
longitudinali), “L’Energia Elettrica, 1934, Vol. 11. 

6. “Ricerche sperimentali sogli sfioratori longitudinali, “L’ Energia 
Elettrica, 1938, Vol. 9. 

7. “Sicherung der Werkkanale gegen eindringendes Hochwasser,” 
Wasserkraft und Wasserwirtschaft, 1941, Vol. 1, p. 12. 

8. “Beitrag zur Abfuhrberechnung von Streichwehre,” Wasserkraft, 1925, 
pp 41, 59, and 99. 

9. “Contribution 4 l'étude des courants liquides,” Zurich, 1933, : 

“Expériences sur le mouvement permanent de l’eau Gane les canaux dé 

couverts, avec apport ou prélevement de long du courant,” Bullen 

technique de la Suisse Romande, 1937. 
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W. H. R. Nimmo!1 applied momentum equation to compute the water level 
ariations in front of the lateral spillway crest. 


2, Nature of Flow Along the Lateral Spillway 


Fundamentally, the flow along the crest of a lateral spillway may be 
treated as a flow division where the divided water, lost over the crest of the 
spillway, does not affect the energy head. The same conclusion may be drawn 
when the flow in the spillway section is regarded as an imaginary expansion. 
Thus, the Bernouilli equation combined with the continuity equation is applic- 
able in computing the lateral spillway action. 

The variation of the water depth, Ay, along a spillway section of the length, 
4x, may be computed easily by the cut-and-try method, using Eq. 1, 


oe (1) 
al 2 
Ay + (S, = S)Ax = Ean © 
f @ -C h¥#Ax 
in which eis rw (2) 


In Eqs. 1 and 2, S and & are bottom and friction slopes, respectively; Vy and 
Vp are the corresponding mean velocities in the upstream and in the down- 
stream end of the section; g is the gravitational acceleration; Q) represents 
the rate of flow in the upstream end of the section; hm is the average head on 
the spillway crest in the section; Ag is the flow area in the downstream end 
of the section; C is the discharge coefficient of the spillway; and Cy is the 
coefficient of the velocity distribution (Coriolis). 

. Eq. 1, friction term S; equalized to the bottom slope S, Cy assumed to be 
unity, and Vo replaced by its value (Eq. 2), can also be derivated from the 
momentum loss in the section. The momentum in the upstream and in the 


downstream ends of the section, Ax, are: 


: My = 4 QV, and 
| Me weep ee ( Questa eh) Vv 
; Dee erly m 2 


in which w is the unit weight of water. Replacing V2 by Vi- AV, the momen- 
tum difference in both section ends, AM 4 = Mj-Mg, can be expressed as fol- 


lows: 


AM, = FL GAY iat wAx( = av) (3) 


‘0 this momentum difference, the momentum lost by the discharge over the 
illway crest: 


1. “Side Spillway for Regulating Diversion Canals,” ASCE Transactions, 
Vol. No. 92, 1928, p. 1561. 
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4 
"6 ni Ax(v, - AY (4) 


has to be added. Thus the total momentum change, AM, in a lateral spillway 
section is expressed as follows: 


Y AV (5) 
M = Z(Q,Av - ¢ We Ax) 


This corrected momentum loss along a section of a lateral spillway, modified 
and equated with the decelerating force in the section in the same way as done 
by Julian Hinds,12 yields the equation: 


¥, 
Q3 % toa C he Ax 6) 
Ay = ze cere — rT 


Eq. 6 is fully identical with Eq. 1 as can be easily shown by replacing AV wit 
3/ 

V1-Vg and C — Ax with Q;-Q9, when in Eq. 1 the friction factor is omitted 

and Cy is unity. 


3. Diversity of Flow Surface Profiles 


The water surface profiles upstream, in front, and downstream from a 
spillway depend on the flow stage, on the relative height of the spillway crest 
on the velocity distribution in the cross-section, and on the location of the 
spillway in the canal. In case of subcritical flow in a canal of mild slope, the 
water depth at the downstream end of the spillway is determined by the con- 
veyance depth in the canal or by the backwater from a control structure or 
from a close forebay. The depth at the upstream end of the spillway is con- 
trolled by the constant energy head along the spillway. Consequently, the 
surface profile upstream from the spillway has to be a straight or either a 
drop-down or a backwater curve, depending on the necessity to maintain, to 
absorb, or to preserve the energy in the approaching flow so that it would 
match the energy head downstream from the spillway. 

Following characteristic variations of the surface profiles may be dis- 
tinguished: - 
a. The lateral spillway located close to the intake gate, Fig. la. Its purpose 

is to spill off the excessive rate of the diversion flow during the high head 

water. The energy loss of the flow, passing through the gate opening, » 
should absorb the excess energy head between the flood level in the reser: 
voir and the maximum permissible conveyance energy head in the canal. 
The water level between the gate and the spillway is affected by the un- 
equal velocity distribution behind the submerged gate opening and in the . 
spillway section. 


V 


b. The spillway is located in the middle of a long canal to skim off the exc 
sive surface or the side-stream inflow, Fig. 1b. The water surface 


12, “Side Channel Spillway,” ASCE Transactions, Vol. 87, 1926, p. 894. 
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downstream from the spillway is controlled by the maximum permissible 

conveyance depth as in the preceding case. As the energy head in the up- 

stream flow is higher than downstream, the approaching flow follows a 

dropdown curve in order to destroy the excessive energy by friction, Of 

course, enough canal length on the upstream side is required to develop 

the necessary down-drop; otherwise, the spiliway would not be as effective 

as intended. 

. The spillway is located at the end of a canal close to the forebay, for the 

purpose of discharging the full inflow at abrupt shutdown of the plant, 

Fig. lc. (The forebay is assumed to be large enough to be able to dampen 

| the closure surge). In this case, the water level at the end of the lateral 

_ spillway is higher than the conveyance head in the conduit. To match the 

higher energy head, the water level upstream from the spillway builds up 

a backwater profile to obtain the necessary rise of energy head. 

The flow in a canal of subcritical slope may pass into the supercritical 

stage at the beginning of the spillway and jump up again to the subcritical 

stage downstream, Fig. 1d. Such undesired flow may occur: 

1. When the spillway crest is below the critical depth.(Prof. A. 
Schocklitsch!3), or 

2. When the ratio of the velocity upstream from the lateral spillway to the 
critical velocity exceeds 0.75 (M. Schmidt14), Under such flow condi- 
tions, the profile of the water surface and the performance of the spill- 
way cannot be determined without model tests. The jump may occur 
already in front of the crest. Further, the spillway capacity cannot be 
increased much by lengthening its crest because the hydraulic jump is 
shifted downstream when the spillway is made longer. 

. Ina canal of supercritical slope, the water surface drops in front of the 

spillway crest, Fig. 1c. With increase of the velocity, the energy level in 

the spillway section cannot be held parallel to the bottom because of a 

strong increase of friction loss as observed by Gentilini.6 


oT 


, Cross-Sectional Velocity Distribution 


As shown in paragraph A-2, the water surface in front of a lateral spillway 
an be computed by Bernouilli equation, imagining that the flow section is 
idened correspondingly to the flow division over the spillway crest. 

- Actually, only the flow section above the spillway crest can expand freely. 
onsequently, an unequal velocity distribution will ensue in the spillway sec- 
on. This condition may be aggravated by the unequal velocity distribution 
shind the submerged gate opening. 

For the computation of the water level, a velocity coefficient, Cy, is intro- 
iced to take care of the depression of the water level caused by the unequal 
-oss-sectional velocity distribution. The magnitude of such velocity coef- 
cients could be established by model tests or by prototype measurements. 
ortunately, two papers concerning the surface depression in the lateral 
jillway flow have been recently published by Dr. Engs. M. Schmidt14 and 


. “Handbuch des Wasserbaues,” Julius Springer, Vienna, 1950, Vol. 1, 


p. 136. 
. “Die Berechnung von Streichwehren,” Die Wasserwirtschaft, Vol. 1, 


1955, p. 96. 
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Zschiersche:15 Cy-values of the magnitude up to 1.3 at the beginning of the 
spillway have been established. Even higher values have been found at the 
end of the spillway crest. 

According to Schmidt, the lateral spillway flow is not controlled complete 
ly by the Bernouilli equation, To adjust the Bernouilli equation for the later 
spillway, Schmidt introduces a ratio, Cy, for the variation of the depth to the 
variation of the kinetic energy along the spillway: 


LEE OS ag ES (7 
1.i 5, - 1-1 58-1 Se : 


Y1and yg are the respective actual depths in up- and down stream ends of 
the spillway; Vj and V9 are the corresponding mean velocities; and LS is 
the friction loss. In Fig, 2a, the values of Cy for two different lateral spill- 
way layouts, obtained from experimental study by Schmidt, are plotted as 

Dn 
z+h 
The curve-(Q9>0) represents the Cy-values for the regular lateral spillways 
behind intake structure as described in case a in paragraph A 3. The curve- 
(Q9=0) is applicable for the lateral spillways located at the ends of the con- 
duits (case c, paragraph A 3). Zschiersche, conducting similar tests, estab- 
lished correction coefficients for empirical use. 

As described previously, the flow over the lateral spillway crest is a flow 
division. Fundamentally, no appreciable change in the energy level can occur 
in a flow division. The higher friction loss due to the increased non- 
uniformity of the cross-sectional velocity distribution is taken care of by 
multiplying the conventional friction loss, LSs, by the coefficient Cy. As men 
tioned above, the friction loss can be compensated by the corresponding bot- 
tom drop. Consequently, the simplified Bernouilli equation: 


functions of the ratio, , in which z is the height of the spillway crest. 


_— - Mere (8) 


is applicable for the computations of the lateral spillway flow. In Eq. 8, Cy 
is 1,1 times the respective Cy-value obtained from the graph in Fig. 2a and 
Hg is the constant energy head in the spillway section. Actually, the velocity 
coefficient Cy varies along the spillway crest (Schmidt). As there is no ex- 
perimental data about its variation available, the value of Cy is assumed to 
be constant along the crest of a particular spillway. For conveniency, the 

h 


Cy-values are plotted as functions of the ratio, . in Fig. 2b. 


m 
Pell ee 
m 


15, “Die Berechnung von Streichwehren auf Grund von Modellversuchen 
geraden und schrag gestellten Streichwehren, “Veroffentlichung der 
Forschungsanstalt fur Schiffahrt, Wasser- und Grundbau No 2, Akademie 
Verlag, Berlin, 1954. 
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According to the experimental studies of the action of lateral spillways by 
ntilini,® the discharge coefficients established by Bazin for regular sharp- 


rested weirs may be also used for lateral spillway crests, if reduced 7 per 
cent: 


Q=cLwH2 in which (9) 
0.0 hs 
C = 0.93(3.248 + sone gees =o (9a) 
Y2 
fee nits: 3,7, eoiann ey Seto yi Forcheimer‘ and Schmidti4 
s HS 


oroposed to use for lateral spillways of any crest shape the corresponding 
lischarge coefficients of the regular weirs reduced by 5 per cent. 


B. Lateral Spillway at Conduit Intake 


|. Performance of Simple Lateral Spillway 


The writer investigates the performance of the most common lateral spill- 
way located behind the intake structure (Fig. 3a). The purpose of this spill- 
vay layout is to protect the canal banks against over flooding during the water 
evel rise in the reservoir. The spillway crest is placed so high that no water 
Ss wasted under normal flow conditions. 

In a uniform spillway section, the water depths can be computed using the 
3ernouilli equation: 


2 
¥ + C, 3, +(Sp- S)Ax = Hy to 


n which the mean velocity, V, is reduced in the downstream direction be- 
ause of the loss of water over the spillway. Replacing V with its value 


y and equating the friction slope S_ to the bottom slope S, Eq. 10 for a uni- 


orm rectangular section is modified to: 


o- EP ay) 


id. 11 is the most important equation in designing lateral spillways in flumes 
f uniform rectangular cross-sections. (Analogous equations for other uni- 
orm canal shapes may be derived in the same way). For an established 
lume width b, velocity coefficient Cy, and energy head Hg at the flood stage, 
1e rate of flow as a function of the depth can be computed. Plotting this data 
raphically in a y,Q-coordinate system yields the (y,Q)-curve. 

_ As an example, the action of a lateral spillway for a 10-ft wide concrete 
ume with a slope of 0.00126 and with a roughness coefficient n (Manning) 
qual to 0.015 is investigated (Fig. 3a). The normal conveyance depth, yn, is 


; 
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6.0 ft for Qn equal to 500 cu ft sec. At first, the velocity coefficient, Cy. is 
assumed to be equal to one. The maximum permissible depth, y¢, in the 
flume is 7.00 ft at 611 cu ft sec discharge and at an energy head of 8.19 ft. 
The elevation of the spillway crest, z, is assumed to be 3 in above the norma 
water level to prevent the water from splashing over it at normal flow. The 
intersection point of the crest line with the (y,Q)-curve in Fig. 3b represents 
the maximum theoretical diversion rate of 697 cu ft sec which the spillway 
may control at 0.75 ft maximum spillway head. However, the length of such 
spillway, when the water surface is level with the spillway crest at its begin- 
ning would be very long and the spillway itself inefficient. The reason for 
that is that the rise of the water level would not be initiated by the discharge : 
action of the spillway, but only by the friction effect in its conduit. To start 
up the rise of the water level in the spillway section, a certain minimum heac 
on the crest at its beginning is required. 

As indicated by the graph in Fig. 2b, the velocity coefficient, Cy, in a spill 
way section may assume values up to 1.33. Therefore, neglecting the effect 
of the uneven velocity distribution in the spillway section may result in inade- 
quate performance of the spillway. At first approximation, the magnitude of 
Cy may be determined using the depths obtained from (y,Q)-curve with uni- 
form velocity distribution. Using the maximum crest head, hg=V2-z, corres- 
ponding to the maximum conveyance flow rate in the canal, and the assumed 
minimum head, hj, the average crest head, hm, is determined. The approxi- 


mate value of Cy obtained as a function of —— from the graph of Fig. 2b 
m 

is somewhat smaller than the final Cy-value. The final Cy can be obtained 

repeating the same procedure with the corrected maximum depth. To be con 

sistent in using the velocity coefficients, the normal water levels in spillway 

sections are computed with a Cy-value equal to 1.05 in all spillway examples 

in this paper. 

In Fig. 4, the upper branch of the same (y,Q)-curve as in Fig. 3b is plottec 
in a larger scale as curve (1). The section of a heavy line below the thinner 
(y,Q)-curve represents the corrected water surface profile on the spillway 
crest with Cy equal to 1.17. As indicated in the Fig. 4, the maximum diver-_ 
sion flow rate at the flood water stage in the 10-feet-wide spillway section is 
only 634 cu ft sec at a minimum depth of 6.42 ft (0.25 ft above the spillway 
crest), The maximum depth at the end of the spillway for the maximum pers 
missible flume flow of 611 cu ft sec at the flood stage is 6.67 ft. The ne 
mum head onthe spillway is 6.67 - 6.17 = 0.50 ft. 

As reported in paragraph A-3d, the approach velocity upstream from thal 
spillway should not exceed 0.75 times the critical velocity, i.e. the theoreti- 
cal depth (at Cy = 1.0) at the beginning of the spillway should not be less than 
1 1/3 times the critical depth. Otherwise, the flow may pass into the super- 
critical stage in spite of a subcritical bottom slope. As indicated by the in-~ 
tersection of (1 1/3 yer)-curve with the (y,Q)-curve, it is quite close to the 
maximum diversion rate but still on the safe side. 

The required length of the side spillway crest, L, is obtained by the au 
mation of the discharges over the sections of the crest, ‘ 

y 


=3 ae (12) 
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: 1 
he conjucated values Q = Q ~ Qo and ee 5} in the crest sections can 


fe easily computed by using Eq. 11 or scaled from the (y,Q)-curve in the 
lange of the spillway. 


. Efficiency of Lateral Spillway 


The best criterion of the efficiency of a side spillway layout is the maxi- 
num rise of the headwater obtainable during the flood stage without encroach- 
nent on the minimum freebord of the canal. 

In general, at a certain gate or orifice opening, the drop of the water level 
a an intake structure is about proportional to the square of the rate of the in- 
ow. That means, the greater the spilling capacity of a lateral spillway, the 
igher is the possible rise of the headwater level during flood. 

Using the discharge capacity of the lateral spillway as a criterion of its 
fficiency, the performance of the simple spillway layout, described in the 
receding chapter, is very poor. As indicated by the (y,Q)-curve (1) in Fig. 4, 
1€ gradient of the depth variation in front of this spillway is rather steep. 
sesides that, the depth at the end of the spillway section, yo, is only 0.25 ft 
reater than the depth, yj, at its beginning. Because of this small change of 
epth and of the sharp downward curved section of (y,Q)-curve in the range of 
lis spillway, the maximum diversion rate cannot exceed much the permis- 
ible conveyance flow of the conduit. 

To improve this condition, the velocities in the spillway section are re- 
uced: 

1. By widening the spillway section to 15 ft or 

2. By deepening it by 3 ft. 

1 this way, the range of the spillways will be shifted into the flatter sections 
f the (y,Q)-curves as shown on curves (2) and (3) in Fig. 4. The energy 
eads, H,, during the flood flow are increased from 8.19 ft to 8.25 ft for the 
idened and to 11.26 ft for the deepened sections. As marked in Fig. 4 and 

1 Table 1, the safe maximum inflow is increased from 634 cu ft sec in origi- 
al 10 ft wide spillway section to 799 and to 805 cu ft sec for the widened and 
wr the deepened spillway sections, respectively. However, the maximum 
jillway head, hg, is increased from 0.50 ft only to 0.74 and to 0.79 ft re- 
ectively. Therefore, the lengths of the lateral spillways would be longer: 
90 ft for the 15-ft-wide and 184 for the 3-ft-deepened sections. 


Improvement of Performance of Lateral Spillways 


As shown, widening or deepening of the spillway section would increase 
e spilling capacity of the lateral spillway. However, the spillways would 
en be very long because of the small heads on their crests. In order to re- 
ice the length of the spillway, an additional orifice section may be placed 
tween the spillway and the flume as shown in Fig. 5. This arrangement, 
sing the same flow conditions in the flume as before, increases the perman- 
tt head loss during the normal flow, but cuts down the length of the spillway 
indicated in Fig. 6 and in Table 1. Besides that, the length of the spillway 
n be shortened when the minimum head at the beginning of the crest is in- 
eased from the minimum 3 as used in all described layouts. However, the 
aximum diversion rate and therewith the possible headwater rise would 


en be reduced, too. 
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To demonstrate thé efficiency of the orifice spillway, 3 alternate layouts 

are investigated: 

1) An orifice of 53 sq ft net area is placed between the flume and spillwa 
sections, both 10 ft wide. This orifice increases the energy head in th 
spillway section to 9.0 ft at the maximum flood flow when the water 
level in the flume is permitted to rise one foot; 

2) An orifice of 43.4 sq ft net area is used to raise the energy head to 10. 
ft. The flume width is 10 ft as in (1); and 

3) The section of the spillway is widened to 15 ft. The flood energy head 
and the net orifice area are the same as in (2). : 

The water level drop, yo-yf, during the flood stage behind the orifice can 

be computed by the following equation: 


2 2 
ee = ee V, 
N Pies INE se (Q/C bs CyVo (13) 


modified from the standard orifice equation.16 In Eq. 13, CgAg represents 
the net orifice area, assumed to be constant for normal as well as for flood 
flow; yo and yg are the respective water depths in front and behind the orific 
and V, is the mean approach velocity. 

In comparison to the permissible flood diversion of 634 cu ft sec in a sim 
ple spillway layout of 10 ft width and of 8.19 ft energy head, the diversion rat 
at the flood stage, when orifices are used, is increased to 710, 825, and to 
1155 cu ft sec, respectively, as shown in Fig. 6 and in Table 1. The most in 
portant advantage of an orifice spillway is the increased head on its crest. 
The maximum head at the end of the spillway is increased from 0.50 ft in the 
original 10-ft-wide spillway section to 1.21, 1.63, and 1.70 ft in the orifice © 
spillways (1), (2), and (3), respectively. As a result, the respective lengths 
of the lateral spillways are 55, 98, and 334 ft. 

As demonstrated by the 15-ft-wide orifice lateral spillway, the flood di-. 
version rate could be increased considerably when the spillway section is en 
larged. However, it is not advisable to use this extreme diversion rate be- 
cause of too long and partly inefficient spillway crest. It is expedient to 
increase the minimum head on the crest and compensate the difference in th 
headwater rise by reduction of the net orifice area behind the spillway. In — 
this way, the necessary crest length can be reduced. : 

As is recognizable from Fig. 6, the surface depression in the front of the 
spillway crests of the orifice spillways; illustrated by the heavy lined curve: 
is much more distinct than in the simple spillway layouts. However, this is 
only accidentally so, because of very unfavorable head-depth ratios in front 
of the spillway crests which effect the most non-uniform velocity distribu- 4 
tions with Cy-values 1.3 and over (Fig. 2b). “ 

The stability of the flow in the orifice side spillways is good as shown by 
(1 1/3 y¢)-curves in Fig. 6, which intersect the (y,Q)-curves at Q-values wi 
above the maximum diversion rates. ; 

A disadvantage to the lateral spillways with an additional orifice is the 
creased permanent loss at normal flow. > 


16. “Handbook of Hydraulics” by H. W. King, 1955, p. 3-10. 
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. Application to Design Problem 


In actual design, the headwater rise, AH, in a reservoir during the flood 
e is determined by the characteristics of the drainage and storage areas 

d by the type and capacity of the discharge structures, Likewise, the free- 

ord of the canal is established with concern of the economics. 

As shown in previous chapters, a higher headwater rise requires a higher 
illing capacity of the lateral spillway. Besides that, economics of the 
ructure and the value of the permanent head loss in intake and in inserted 

+ ia have to be considered, Further, the flow in the intake and in the later- 
1 spillway sections should not pass into the supercritical stage. 

: In preliminary studies, it is useful to employ a simplified equation for es- 
ablishing the required net area of the intake gate opening, CjAj; 


2 2 
sorry Sib Be athe 7 
Mi \V Be(AH-ay,) ) 


1 which the constant orifice coefficient, C;, is assumed to take care of all 
asses in the intake at any flow rates. Qp andAy , =Y1-Yn, are conjugated 
alues established by Eq. 11. 

As shown by Eq. 14, it is theoretically possible for any rate of floodwater 
itake, Qp, and headwater rise, AH, to establish a net gate opening. How- 
ver, there are physical and economical conditions which limit the size of 
1e gate opening: 

1) At certain Qp-values, CjA; cannot be reduced beyond a minimum area, 

otherwise the flow behind the gate may pass into the supercritical stage, 

2) At smaller gate openings, the permanent head loss during the normal 

flow may be prohibitively high; and 

3) As demonstrated by (y,Q)-curves in Fig. 4, the required length of the 

spillway crest may be too long to be economical. ni 

Further, Eq. 14 indicates that an increase of the water level rise, A yj, at 
1e beginning of the spillway helps to reduce the required maximum diversion 
ow at the same net gate opening. This way, the crest length of the spillway 
in be shortened. 

- Considering the foregoing statements, it is obvious that trials are neces- 
ry in designing economical and safe-acting lateral spillways. Being famili- 
- with the characteristics of the performances of different lateral spillway 
youts, explained in this paper, the designer will be able to cut down the num- 
sr of trials. 


CONCLUSIONS 


In general, the lateral spillway itself is not capable of destroying any of 
e excess energy head upstream from it. The purpose of a side spillway in 
‘otecting an open canal against overflooding is merely to spill off the sur- 
us of the flow. The excessive energy head has to be absorbed upstream 
om the lateral spillway by the friction, intake and gate losses, or under 
rcumstances even by a controlled hydraulic jump. Ina special case, when 
€ spillway is located in front of the power intake to discharge the full or 
rt flow at the shut-down operation of the turbines, no energy head is 
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destroyed, On the-contrary, the necessary energy head, to raise the water 

level, is built up by reduced friction loss in the backwater. 

In respect to the computation and improvement of the lateral spillway, 
placed behind the diversion intake, the following conclusions can be drawn: ~ 
1) The performance of a lateral spillway can be computed using the | 

Bernouilli and continuity equations. Because of the non-uniform velocity — 

distribution in the lateral spillway section, the velocity coefficient, Cy, 

has to be considered in computations. 

2) In uniform sections of the lateral spillway, the variation of the water level 

profile and the maximum spilling capacity, characteristic for each spill- 

way layout, can be computed by Eq. 11. 

It is useful to plot the (y,Q)-curves in designing lateral spillways. The 

relative location of the maximum diversion flow in the curve gives ad- 

vance information about the spilling capacity and about the safety against 
the possible change of flow stage in the spillway section. 

As demonstrated by the performances of different spillway layouts, the 

simple spillway with comparatively high velocities has a small spilling 

capacity. Consequently, the possible headwater rise is very limited. 

Enlargement of the spillway cross-section helps to increase the spilling 

capacity and therewith the capacity to raise the flood level of the head- 

water. However, the spillway would be considerably long because of the 
small head on its crest. 

6) An improvement can be obtained when an additional orifice is placed be- 
tween the spillway and the canal. Although, the permanent loss at normal 
flow conditions is increased by the orifice loss, the length of the spillway 
is shortened considerably because of the higher head on its crest. 


3 


~~ 


4 


— 


) 
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APPENDIX. NOTATION, 


The following symbols, adopted for use in the paper, conform essentially 
with Letter Symbols for Hydraulics (ASA Z10.2 - 1942), 


Ay = area of the wetted cross-section at the end of a section, in sq. ft; 

b = width of the section, in ft; 

Cc = coefficient of spillway discharge; 

Cy = correction coefficient of Bernouilli equation, by Schmidt; 

Cy = coefficient of velocity distribution; E 


CoAg =net orifice area, in sq ft; a 


CjAj = net area of the opening of the intake gate, in sq ft; 


g = acceleration of gravity, in ft per sec per sec; ; 
hy = head on the crest at the beginning of a section, in ft; ; 
ho = head on the crest at the end of a section, in ft; 4 
hm =‘ = average head on a spillway section, in ft; 4 
Hp = energy head of the headwater at flood flow, in ft; 4 


Hy = = energy head of the headwater at normal flow, in ft; 
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= energy head in the spillway section at flood flow, in ft; 


dé tn 


= rise of the headwater at flood flow, in ft; 


rr & 
is) 


= length of spillway crest, in ft; 


= momentum in the upstream end of a section, in sec pound; 


oP 


= momentum in the downstream end of a section, in sec pound; 


AM, = momentum difference between the section ends, in sec pound; 
AM = total momentum change in a spillway section, in sec pound; 

2 = coefficient of roughness in the Manning formula; 

a) = rate of flow, in general, in cu ft per sec; 

ay = rate of flow in the upstream end of a section, in cu ft per sec; 
29 = rate of flow in the downstream end of a sect., in cu ft per sec; 
aN = rate of diversion at normal flow, in cu ft per sec; 

QF = rate of diversion at maximum flood flow, in cu ft per sec; 

AQ = discharge across a section of the spillway crest, AQ = Qi - Q9; 
Sf = slope of the energy head due to the friction loss; 

Ser = critical slope of the canal; 

5 = slope of the canal bottom; 

Jo = mean approach velocity in front of an orifice, in ft per sec; 

4 = mean velocity in the upstream end of a section, in ft per sec; 
9 = mean velocity in the downstream end of a sect., in ft per sec; 
\V = change of mean velocity ina section, V = V,-Vg9, in ft per sec; 
v = unit weight of water, in pounds per cu ft; 

\x = length of a section of the spillway, in ft; 

n = normal conveyance depth of flume, with Cy = 1, in ft; 

£ = flood conveyance depth of flume, with Cy = 1, in ft; 

er = critical depth of flow, in ft; 

y = variation of water level, in ft; 


depth of water in the spillway section, with C,>1, in ft; 


= depth of water in front of an orifice, with Cy>1, in ft; 


° 

1 = depth of water at the beginning of the spillway, Cy>1, in ft; 

3 = depth of water at the end of the spillway, with Cy>1, in ft; 

v1 = rise of flood water level at the beginning of spillway, Cy-1, in ft; and 


elevation of the crest above the bottom, in ft. 
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THE PROBLEM OF RESERVOIR CAPACITY FOR LONG-TERM STORAGE 


A. Fathy* and Aly S. Shukry,** M. ASCE 
(Proc. Paper 1082) 


SYNOPSIS 


Until very recently there was no theoretical basis for the determination of 
the capacity needed in a long-term reservoir to guarantee a given draft. 
Lately, a solution of the problem has been given by Hurst on the assumption 
that hydrological occurrences may be treated as random events the sequence 
of which is free of any periodic trend. Hurst has shown that, subject to that 
condition, the main factor by which the capacity is governed is the standard 
deviation of a single observation. 

Considering that in hydrological phenomena there is usually some tendency 
towards grouping of high or low years, the authors have approached the prob- 
lem from a different angle. In the new treatment, deviations in the arithmetic 
mean for groups of observations are utilised instead of the deviation of one 
observation. In this way, the influence of any peculiarity of sequence in the 
phenomenon under consideration will obviously be covered. 


INTRODUCTION 


Ideal Form of the Problem 


In the design of over-year storage schemes for irrigation, power genera- 
tion or town water-supply an important consideration is that of reservoir ca- 
pacity. In its simplest form, the problem consists in the determination of the 
capacity required to secure a constant draft equal to the mean supply. The 
fact, however, that it is usually impossible to make an accurate forecast of 
the supply for a long time in advance renders the notion of a draft exactly 
equal to the mean supply impracticable. Nevertheless, that notion constitutes 
the logical starting point for a theoretical study of the problem, A general 
idea of the fundamental principles involved will be gained from the following 


2xample: 
————E——— ; 
Note: Discussion open until March 1, 1957. Paper 1082 is part of the copyrighted 
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Suppose we are given a set of observations of variable magnitude repre- 
senting the annual discharge of a river at a certain site, and it is required to 
find the capacity needed in a balancing reservoir with the aid of which we can 
secure a constant annual draft equal to the mean of these observations. 

It is easy to see that the capacity required equals the range between the 
maximum and the minimum values of the cumulated departure of individual 
observations from the mean. This is illustrated in the following table: 


Table I 


Remarks 


wow onrsaw Ff WwW WYP H 


In this case, the storage needed will be the range between 24 and -12, 
which is 36. 

Another way in which this result may be arrived at is to tabulate the 
cumulated draft against the cumulated supply. The difference between the 
two at the end of each year will be the reservoir content at that instant. This 


is shown in the next table. 
A peculiarity in Table II is that at certain times the reservoir ‘5 


< 


contents become negative. A constant draft equal to the mean cannot, there-— 
fore, be obtained throughout the whole of the working period unless we borrow 
a volume at least equal to the maximum negative storage. Supposing a ‘a 
volume of 12 to be borrowed at the start, the adjusted reservoir contents will 
be as shown in the last column. Obviously, the residue in the reservoir at 
end of the operation must correspond to the borrowed volume since the t 
outflow equals the total input. 


Graphical Representation 
The graphic plot of either the cumulated supply or the cumulated draft is. 
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Adjusted 
pou nen 
0) 


Table II 


oo ON A ee AN hd 


H 
to) 


called a “Mass-Curve.” The mass-curves for the case assumed above are 
shown in Fig. 1. The curve Oabc represents the cumulated supply. The 
Straight line Oc represents the cumulated draft from an empty reservoir, 
while line O'c' represents the cumulated draft with an initial storage of 12. 
The following points are useful to note: 


1) The slope of a mass-curve at any point represents the rate of flow. 

2) For a constant rate of flow the mass-curve is a straight line. 

3) The reservoir content at any instant is given by the respective intercept 
between the mass-input and the mass-outflow curves. 

_ 4) The range is given by the intercept between two lines drawn parallel to 
the output line through the highest and lowest points on the input curve in re- 
lation to that line. 

_ 5) The mean up to any point on a mass-curve is given by the slope of a 
straight line drawn from the origin to that point. Similarly, the mean for any 
inter mediate period is given by the slope of a straight line joining the initial 

id the terminal points on the mass-curve, 


The intermediate period determining the range (from a to b in Fig. 1) will 
re be referred to as the “Critical Period.” In the above example, the 
ritical period is a period of deficit. According to its definition, the range 
y also correspond to a period of excess. That depends on the order in 
ch the points of maximum and minimum cumulated departure are en- 
tered. 
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Denoting the length of the critical period by Ne, the respective mean by 
Mc, the overall mean by M,) and the range by R, it is clear that: 


R=N, |u, - u,| (1) 


Variation of Range with Number 


If we take a long series of observations of non-periodic character and 
compute the ranges for sets of different length, we shall generally find that 
the range tends to increase progressively with the number of observations 
taken, Thus, under the conditions assumed above, the capacity needed in a 
long-term reservoir would be indeterminate unless the length of the working» 
period were fixed. A theoretical study of that problem has recently been 
made by Dr. H. E. Hurst in connection with storage schemes in the Nile 
Basin.(1) He has succeeded in establishing a relation between the mean 
range (R), the number of observations (N) and the standard deviation (0), for 
random events of normal frequency distribution thus: 


R = 1.25 o/N 
Noting that (2) 
2 
See (a 
N 


where v = deviation of a single observation from the absolute mean, the above 


expression is equivalent to: 
= 1.25VSv° (3) 


showing that the range increases in the same manner as the accumulated er- 
ror on a line of levelling. 

Finding that the range for natural phenomena did not, in general, confor mi 
to that law, Hurst assumed a similar relation of the form: 


R=ac nt (4) 


where the coefficient a and the index k had to be determined empirically. 
If the best values of a and k were extracted for each phenomenon inde- 
pendently and without, in any way, predisposing the result to be expected, it — 
would be found that the range of variation therein from one phenomenon to E 
another is somewhat large. The result of an analysis carried out by the * 
authors on rainfall observations at seven stations, each with a record of not : 
less than 100 years duration, is given in Table Il. * 
To gain some idea of the influence of this variation on the range, when re- 
lation (4) is applied, let us take N = 100 as a standard for purposes of com-— 
ee For Madras R/o = 9.5. For Boston, R/o = 25.9. 
In dealing with this problem, however, Hurst supposed that the coefficient 
a and the index k in (4) should be such as to satisfy the condition R/o =1 
when N = 2. On this supposition (4) was reduced to the form: 
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Fig. 2 — MAXIMUM DEFICIT 
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Table III 
Station a k 
Copenhagen Feil 0.35 0.79 
Charleston (U.S.A.) 0.54 0.94 
Boston (U.S.A.) 0.13 1315 
Philadelphia (Wie Sista) BRR HOPE hoy 
New York (U.S.A. ) attul 0.74 
Madras (India) I.82 0.36 
Edinburgh (Scotland) 1.82 0.43 
R 
k 
mewn canes (5) 
oy 


Thus the coefficient (a) was tied up with the index (k). By logarithmic plottin 
of R/o against N, in which all lines were made to pass through the point (0, 
log 2), it was found that variation in k was confined within a reasonable range 
for all sorts of natural phenomena such as river levels and discharges, rain- 
fall, temperature and pressure, annual growth of tree rings, mud deposits in 
lakes, etc. For that index an average value of 0.72 was adopted. It followed 
that the coefficient (a) should have the value 0.61, whence relation (4) acquire 
the final form: 


R = 0.61 o N°°/@ (6) 


In the authors’ opinion, that procedure would be justifiable only if the rela 
tion were in the nature of an exact law. Considering that the condition R/o = 
when N = 2 is far from being satisfied in relation (2), for which there is a 
theoretical foundation, it would appear that the consistency apparently gained 
through that extrapolation is the product of the procedure itself and not of an 
uniformity in the character of the various phenomena examined. 


Case of a Draft Less than the Mean 


The question of the storage needed to guarantee a draft less than the meal 
which is the more important from the practical point of view, was considere 
by Hurst on the supposition of the existence of a relationship between that 
storage and the normal range. (1) His treatment of this subject may be 
summed up as follows: 

Let No = overall number, Mo = overall mean, B = draft (less than Mo) anc 
S = storage required to guarantee that draft. The difference between B and 

- Mo was expressed as a fraction of the standard deviation, and S was taken as 
the greatest accumulated deficit, within the period No, with respect to the 
draft B. y 

The notion of the maximum deficit in relation to a given draft is demon- 
strated in Fig. 2. The range (for a draft = Mo) is given by points a andc. 
For the draft B, as represented by line Oe, the maximum deficit is given by 
points a andb, The corresponding storage S must obviously be smaller than 
Re 
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Through analysis of several sets of random and natural observations, the 
following empirical relations were arrived at: 


S _- 
For random events, R= 0.-9I- 0.89 = (7) 
For natural events, Dera oi= 6.05.) 02 (8) 
R o 


It was admitted that a great deal of research was devoted to finding a the- 
oretical relation between S and R but without success. 

Evidently, the first numerical term in either (7) or (8) should have been 
unity since S should have the same value as R when (M-B) = O. The discrep- 
ancy is due to the fact that S has to be taken as the maximum deficit only 
whereas R may correspond either to a deficit or to an excess. 

It is to be noted, however, that once the notion of a draft exactly equal to 
the mean is relinquished, the situation undergoes a substantial change. With 

a draft equal to the mean, the range generally increases progressively as the 

Tumber increases. This cannot always be the case with a reduced draft. In 
fact, if the draft were reduced to the minimum value of a single observation, 
S would be zero whatever the number. Thus there are no grounds for the 
Supposition of a fixed relation between S and R. 


Practical Considerations 


In Art. 1 above it has been pointed out that the notion of a draft exactly 
€qual to the mean is impracticable in view of the fact that that mean could 
not be accurately known in advance. Another reason is that, even if the mean 
were known, some initial storage would, as a rule, be needed to eliminate 
negative reservoir contents and the amount of that storage could not possibly 
be determined beforehand. Thus, regulation would almost certainly have to 
deviate from the ideal program visualised and, once that deviation took place, 
the theoretical structure built up on that notion would tumble down. 

In addition to these considerations, there are two points which have an im- 
portant bearing on the storage problem from the practical side. The first is 
that the capital outlay for a long-term storage scheme is usually high and the 
Capacity needed should be computed on the most specific lines possible. In 
other words, as much regard as possible should be paid to the peculiarities 
of the case in hand. Thus it would be inadvisable to rely on formulae derived 
through taking averages for a large number of phenomena of different charac- 

er. 

The second point is that a small reduction in the draft below the mean is 
usually found to lead to a relatively huge reduction in the capacity needed. 

us we should always allow for some discount in the draft below the esti- 
Fated mean, apart from the allowance for storage losses, even if the whole 

the supply were needed for consumption. That discount would serve both 
as a factor of safety (in the absence of accurate knowledge of the future mean) 
and as a means of securing a great economy in cost with a little sacrifice in 

onefit. 

- With these considerations in mind, the authors have attacked the problem 

rom a new angle, regard being paid in particular to the case of a draft less 
an the mean. It is to be borne in mind, however, that the problem is one 
t cannot be solved successfully by hard and fast rules and that, whatever 
e line of approach, the procedure cannot be freed entirely from the element 
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of conjecture. The most effective tools in the hand of the designer will be 
sound practical judgment and a thorough understanding of the fundamental 
principles involved. 


Characteristics of Hydrological Phenomena 


General Considerations 


The principal phenomena with which we are concerned in the present study 
are rainfall and river flow. In subjecting any such phenomenon to systematic 
analysis, it is tacitly assumed that the occurrence repeats itself regularly 
and that we can rely on obtaining an average annual supply of a certain order 
for an unlimited number of years. Thus, we must exclude from our discus- 
sion such phenomena as occur only at long intervals or in widely varying de-- 
grees of intensity. It has also to be assumed that the occurrence under con- 
sideration, though variable from year to year, is not subject to progressive 
secular change. 

The first question that presents itself here is how far hydrological phe- 
nomena may be treated as random events. One respect in which a difference 
is most likely is that of the frequency distribution. Most phenomena of that 
kind show skew distributions, positive deviations from the mean being gener- 
ally higher and less frequent than negative deviations. This will be readily 
understood when it is noted that while there may be years of more than 100% 
excess above the mean, there cannot possibly be any year with more than 
100% deficit. 

Another distinctive feature in hydrological phenomena is that high and low 
years often occur in runs or long-range cycles. Although these cycles do not, 
as a rule, take place in symmetrical waves or at regular intervals they are 
likely to affect the storage question considerably. In the case of a purely 
random occurrence, it may be taken as an established fact that the normal 
range generally increases as the square root of the number of observations ~ 
when that number is fairly large. The authors have applied Hurst’s formula 
(2) to observations obtained by throwing five dice a thousand times, the sum ~ 
obtained at each throw being recorded as one reading, Although deviations 
on one side or the other took place, the average range agreed very closely 
with the theoretical value. With natural phenomena, however, the average 
range has been found to increase with number in some cases at a slower rate 
and in others at a higher rate than indicated by either (2) or (6). This is due 
mainly to the peculiarities of sequence in such phenomena. = 

Many attempts have been made to prove the existence of regular perio- _ 
dicities in meteorological or hydrological phenomena but, so far, without suc- 
cess. An example is the theory put forward by Brooks that the levels of Lake 
Victoria fluctuate in unison with the number of sunspots, which is known to _ 
change more or less periodically. That theory has been questioned by Hurst 
and Phillips on several grounds. (2) 

With the object of detecting periodicities in the flow of the River Nile, an. 
analysis of the maximum flood levels on Roda Guage, of which the records 
extend back for 1300 years, was made by Jarvis and discussed by Hurst, 
Shumann and others.(3) No definite conclusions were arrived at. It is rea-_ 
sonable to suppose that some correlation must exist between hydrological 
phenomena and periodic cosmical phenomena. The multiplicity of factors 
involved, however, renders any inquiry in this respect futile. 


a 
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A simple method of revealing cyclic trends in rainfall or river flow re- 
cords is to draw curves representing 3-year means, 5-year means, 7-year 
means and so on, each mean value being plotted against the median year of 
the corresponding period. To obtain smooth mean-curves, the subordinate 
periods should not be taken consecutively, i.e. each beginning at the end of 
the preceding one, but should be shifted forward only one or two years ata 
time. Fig. 3 shows a plot for the rainfall at Boston (U.S.A.). For purposes 
of comparison, a similar plot for the 5-dice observations is also given. 

Incidentally, it may be mentioned that if a definite long-range periodicity 
exist, the operation of taking long-period means progressively as suggested 
above will not obscure that cycle, provided no other cycles are present. 
This fact may be demonstrated as follows: 

Suppose we have a function y = sin(nx) and we take means by small shifts 
for a subordinate period of length 2z. The difference between an individual 
observation and the coincident mean will be: 


sin(nx) - if sin(nx) dx = sin(nx) [2 = HAC (9) 
224 x=Z nz 


Thus, the effect of taking means is simply a distortion of the y-scale with- 
out any alteration in the period length. If superimposed deviations are hap- 
hazard, their effects on the mean will cancel out. 


The Mean 


Before discussing the characteristics of the mean in natural phenomena, 
it may be useful to review these characteristics briefly as regards purely 
random events. The latter are generally subject to the following conditions: 


1) Positive and negative departures from the mean are equally likely. 

2) There is usually a theoretical or absolute value of the mean. In the 
case of the 5-dice, for example, the maximum possible reading is 30 and the 
minimum is 5. The absolute mean, therefore, is 17.5. 

3) The absolute mean is most likely to fall within the range of observed 
values, provided the number of observations is not very small (not less than, 
say, 10 or 12). 

_ The significance of these conditions is perhaps most clearly brought out in 
he case of the repeated measurement of a certain dimension such, for exam- 
le, as the distance between two points. If there is no systematic error in 
he process of measurement, we may expect that beyond a few trials the true 
ength sought will lie within the range of the observed values and not very far 
rom their arithmetic mean. In other words, it will be most improbable that 
ll the observed values will be greater or less than the true value. 

Natural phenomena generally differ from random events in all of these 
‘espects. The absolute mean is always indeterminate and individual values 
nay lie above or below the overall mean adopted for long periods on end. 
the condition of equal probability of positive and negative deviations is also 
ot satisfied in the majority of cases. Under these circumstances, there will 
no guarantee that the mean for a given number of future years will not 
leviate sensibly from the mean derived from past observations, On this ac- 
t, consideration of possible variations in the mean becomes a matter of 


amount importance. 
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A study of fluctuations in the mean annual rainfall at a large number of 
stations situated in different parts of the globe was made long ago by Binnie. 


- Binnie.(4) His main conclusions are: 


1) To obtain a reliable value of the mean, the existing record should ex- 
tend over a period of at least 35 years. This was supposed to yield a result 
correct to within 2% of the mean for any longer period. 

2) The extreme positive deviation ina subordinate-period mean from the 


overall mean is generally greater than the extreme negative deviation. 


3) The range of fluctuation in the mean is generally greater for stations of 


low rainfall than for those of heavy rainfall. 


As regards the first point, it may be mentioned that in calculating the 
subordinate-period means consecutive periods only were taken. The ob- 
served maximum deviations were thereby greatly under-rated owing to the 
consequent reduction in the number of periods tested as the period length was 
increased. To draw the fullest possible information from the available re- 
cords, overlapping periods should have been taken, as suggested in the pre- 
ceding article. 

The authors have made an analysis of rainfall records at ten stations, 
each of not less than 100 years duration, to find out the effect of increasing 
the number of shifts for a given period-length. The result is shown in Fig. 4. 
In the compilation of that figure, the average extreme deviations for the ten 
stations were taken. For the sake of uniformity, the subordinate periods 
were shifted five years at a time and the percentage deviation was referred 
to the mean of the first period (which would be the only one known when an 
attempt is being made to estimate the future mean). This analysis shows 
that much greater deviations than those found by Binnie are possible, even 
with periods exceeding 35 years. 

As regards the second point, the cause of the distinction between positive 
and negative deviations has already been referred to in the preceding article. 

The third conclusion indicates that heavy rainfall is equivalent to a com- 
bined supply from several sources between which there is no correlation. 
This may be demonstrated as follows: 

Let there be two exactly similar sources, each having a mean M anda 
standard deviation ao. The ratio ¢/M may be taken as a measure of relative 
‘deviation for one source. For the two sources together, the resultant mean = 
2M and the resultant standard deviation =/20. Thus, the relative deviation 
for the combination will be 0.707 of the relative deviation for one source. 

The question of the maximum deviation in the arithmetic mean, to which 
reference has been made before, is one of major importance in the present 
study. A clearer view of its significance will be gained if we discuss it first 


from the theoretical point of view. 
For random events, it is known that the probable deviation of the arithme- 


tic mean is: 


rt Ce 


This represents the deviation with respect to which the actual deviation is 
equally likely to be greater or smaller in magnitude. It is also known that 
the probability of a deviation greater than a certain multiple of the P.D. 


3 


: 


1082-12 HY 5 October, 1956 


gatas 
VaEaaen 
WA 
Number of Shifts 


Percentage Extreme Deviahon 


Fig. 4 — PERCENTAGE EXTREME DEVIATION 
IN SUBORDINATE=PERIOD MEANS 


(51) 


eo ee eee 


(N) 
Fig. 5 — DEVIATION-CURVE FOR RANDOM EVENTS 


ASCE FATHY - SHUKRY 1082-13 


decreases very rapidly as the proportion increases. Some values of that 
probability are given in the following list:(5) 


Ratio to P.D. Probability 
See Sate et le le EE a 9 Ta 

I.00 1 
2.44 I/10 
3.06 1/25 
3-45 I/50 
3.82 I/100 
4.85 I/1I000 
5.86 I/10000 


It is quite logical to suppose that probability conditions with respect to the 


: arithmetic mean (regarded as an observation in itself) must be similar to 
those relating to single observations, with the difference that the standard 


deviation for the former = 0/ VN instead of o. Thus we can say that the 
probability of meeting with a deviation greater than, say, 3.45 times the P.D. 
in a set of observations of the mean for some fixed number of individual ob- 
servations = 1/50. Or, to put it the other way, if we have 50 values of the 
mean for a given number N then, on the average, one of these will depart 
from the absolute mean by 3.45 times the P.D. or more, on the understanding 
that the P.D. is that corresponding to the number N. 

Although the deviations represented in Fig. 4 are referred to the mean of 
the first period and not to the absolute mean, the plot shows clearly that na- 
tural phenomena are subject to the same general trends. It is also to be 
borne in mind, however, that variations in the mean for such phenomena are 
most likely to be influenced by peculiarities of sequence as well as by the 
standard deviation. On this account, it is not considered admissible to apply 
the above mentioned rules, without reservation, to any natural occurrence. 
Further reference to this question will be made later. 


The Standard Deviation 


Just as we cannot assign an absolute value to the mean for natural phenom- 
ena, so is the situation as regards the standard deviation. Observation shows 
that in such phenomena the range of fluctuation in the standard deviation is 
generally much wider than the range of fluctuation in the mean. For example, 
the means and standard deviations for the rainfall at Greenwich in three con- 
secutive periods, each of 30 years duration (1841-1930), are: 614, 602, 629 
and 120, 83, 108, respectively. In these figures, the highest value of the mean 
is 4.5% greater than the lowest, while the highest value of the standard devia- 
tion is 45% greater than the lowest. Another striking example is furnished by 
the annual mud deposits in Lake Sake (Crimea). For two consecutive periods, 
each of 100 years, the means and standard deviations are: 10.9, 15.8 and 3.9, 
21.0, respectively. In this case, the higher value of the mean is 45% greater 
than the lower, while the higher value of the standard deviation is 440% great- 
er than the lower. These figures are taken from tables attached to Hurst’s 
paper already referred to,(1) According to Hurst, the average range of vari- 
ation of the mean rainfall, for periods of 40 to 50 years, is 14%, while the 
corresponding range for the standard deviation is 23%. 


I 


] 
. 


1082-14 HY 5 October, 1956 


The Range of Cumulated Departures 


A study of the manner of variation of the range of cumulated departures 
with number will provide us with a guide as to how far natural phenomena 
may be regarded as truly haphazard, As Hurst has established theoretically, 
the mean range varies directly as the square root of the number, in random 
events of normal frequency distribution. The coefficient of proportionality 
depends on the standard deviation but, as far as the general behaviour of the 
range is concerned, it is immaterial on what that coefficient depends so long 
as it is a constant for the phenomenon considered. 

From tests carried out by the authors on the direct relation between range’ 
and number for natural phenomena, it appears that the range generally in- | 
creases with number more rapidly than for random events, the exponent of N 
in some cases being very near the value unity. Owing to the shortness of the 
records, however, no general conclusions can be drawn from these tests. 

In this respect, the authors are inclined to the view that the relation be- 
tween range and number is influenced not only by the standard deviation but 
also by the peculiarities of sequence which are observable in almost all 
natural occurrences. This will be readily understood when it is noted that for 
a perfectly periodic phenomenon the range has a limiting value which is never 
exceeded no matter how great the number. 

Irregular cycles, on the other hand, may have a reverse effect on the 
range, i.e., may cause it to increase with number more rapidly than when 
the distribution of individual observations is absolutely random. On this ac- 
count, it is believed that every phenomenon should be dealt with according to 
its own particular circumstances. 


New Line of Approach to the Problem 


Basis of Proposed Treatment 


In the proposed treatment of the long-term storage problem attention is 
paid in particular to the case of a draft less than the mean, which is the only 
one of practical importance. For purposes of demonstration, however, we 
shall start with the ideal case of a draft exactly equal to the mean. 

Going back to Eq. 1, it will be seen that the range is simply the outcome of 
deviation in the arithmetic mean for some intermediate period from the over- 
all mean for a given number of observations. If the phenomenon is supposed 
to possess an absolute mean (Mg), the overall mean (Mo) will rarely coincide 
with Mg, despite the fact that the latter is the “most probable” value of the — 
mean for any period. Hence we can say that the range is governed by three 
factors, namely: 


1) The length (N,) of the critical period. 
2) The deviation of the critical-period mean (Mc) from Mg. 4 
3) The deviation of Mo from Mg. 


If 6 is the relative deviation of any period mean with reference to the ab- 
solute mean, it follows from (1) that: 


3 
fa x m8, - 3,| a) . 


- 
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For a given value of No, there are many possible combinations between 
Nc, 59 and 6,. Each of these factors will have a probable value of its own. 

It is obvious, however, that we cannot obtain the probable value of R/Mag di- 
rectly by inserting the probable values of the above mentioned factors in (10). 

There may be a means of elaborating that problem on theoretical consider- 
ations but it is not considered worth while to make that attempt because the 
probable range itself is of doubtful practical value. The adoption of the prob- 
abie or the average range would mean that the chances of satisfying or not 
satisfying the prescribed conditions would be about the same. In other words, 
it is the maximum and not the mean range that we should endeavour to ascer- 
tain. 

One way of dealing with that problem is to adopt a standard value for the 
ratio N¢/No and lay down some empirical rule for the choice of suitable 
values of 6, and 6,. This is the essence of the proposed method of tackling 
the problem. In the ideal case visualised, this procedure may appear some- 
what arbitrary but the method is particularly suited to the case of a draft less 
than the mean, as will be shown presently. 

In the application of that method, the essential natural characteristic of the 
phenomenon under consideration to be utilised is the maximum deviation of 
the arithmetic mean for any given number from the absolute mean. That 
characteristic has already been considered for random events in Art. 7. As 
regards hydrological phenomena, it will be assumed that we can construct a 
curve representing the maximum relative deviation that may be expected, 
within reason, to take place in the mean for a given number, with reference 
to a hypothetical absolute mean. 

To make clear that notion, we shall first show how such a curve may be 
constructed for a random event. 

In Art. 7, a scale has been given of the greatest deviations that are most 
likely to be met with in various numbers of observations of the arithmetic 
mean. Theoretically, there is no limit to the relative deviation that may be 
attained. In practice, however, the magnitudes of individual occurrences are 
usually confined between certain physical limits beyond which the mean can- 
not possibly pass. But the probability of attaining any of these limits with 
sets of more than 3 or 4 observations is practically nil since that could hap- 
pen only if all the readings coincided with either the maximum or the mini- 

"mum value. 

Considering that for deviations in the arithmetic mean exceeding a few 
times the P.D. the probability becomes very small, an arbitrary limit may 
be set at some probability which is considered low enough. Supposing we 
take that limiting probability = 1/50, the corresponding value of the most 

likely maximum deviation may be taken = 3,5 times the P.D. On this basis 


we may write: 


2.36 
Max oDe = Me eee 6 
. = (11) 
? Pld. ee 3 OtOy 
anal§] = Si = (12) 


: } a 


The plot of 16 | max, 2gainst N will be called a “Deviation-Curve.” From 

' it i é i the two depends on the ratio o/M 

_ (12) it is clear that the relation between the two dep 4 
which generally varies from one phenomenon to another. Fig. 5 shows a plot 
_of the quantity 6M,/o against N. 

; 
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The next questiom to consider is the relative length of the critical period. 
As already pointed out, the critical period may be either one of deficit, as 
shown at (a), Fig. 6, or one of excess as shown at (b). In the case shown at 
(c) in the same figure, either the excess period Oa or the deficit ac may be 
regarded as the critical period. 

It would be natural to suppose that if deficits and excesses were equally 
likely, the most probable value of Ne would be No/2. This probability is 
most apparent in the last case. Analysis of the 5-dice observations, however, 
shows that the average value of the ratio Nc/No is somewhat less than 1/2 
and that it tends to decrease slowly as No increase. For numbers between 
20 and 200, that average has been found to conform closely to the relation: 


-O.1 

Nosy = 0.67 N. (13) 

This result gains some weight from the fact that the same relation has 
been found to hold for natural as well as for random events. It is evident, 
however, that a highly refined treatment of that point would be superfluous. 
It is, therefore, proposed to adopt for that ratio a standard value of 1/2. 

Finally we come to the question of what values of 69 and 6, we have to 
adopt in relation (10). It stands to reason that if we assign to 6¢ its maximum 
value, as previously defined, we can assume 6, = 0 on the ground that posi- 
tive and negative values of that deviation are equally likely. Strictly speak- 
ing, that argument is not quite correct because Ne is a fairly large fraction 
of No and any deviation in Me from Mg, must to some extent be reflected in 
Mo. In other words, the most probable value of 69 in combination with the 
maximum of 6, would be a small quantity of the same sign as ¢. That, 
again, is a refinement for which there is no need. Relation (10) may now be 
reduced to the form: 


ete ACE (14) 
a 


on the understanding that Nc = No/2 and 4, is the maximum for the period 
Ne. 

Substituting for Nc and 6¢ is (14) their equivalents in terms of No and 0 
we get: 


Pernt fb ov, (15) 


Comparing (15) with (2), we see that the maximum range according to the 
foregoing procedure is about 33% greater than the mean range. This result 
agrees well with actual tests carried out on the 5-dice observations. 

The Deviation-Curve for Hydrological Phenomena ; 

The deviation-curve for a natural phenomenon is obtainable only through ‘ 
direct analysis of the existing observations. It could be readily derived from — 
such a plot as shown in Fig. 3. It is admitted that if the available recordis 
short, the construction of that curve will be difficult. But, in that event, it — 
will be equally difficult to make a reliable estimate of the standard deviation 
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: which is the only alternative to employ. In any case, some judicious adjust- 
_ ments to the observed maxima of the deviations in the means will be neces- 
sary. 

The maximum deviation for a given number of observations may generally 
_ be expressed by a relation of the form: 


Piles ome (16) 


: N 


where the best values of c and m have to be determined for each phenomenon 
independently. If the available record is short, then we may either use the 
theoretical relation for random phenomena (II) or adopt some relation derived 

for another record of long duration of which the general characteristics are 
believed to be similar to those of the phenomenon under consideration. 

As there are only two unknowns in (16), the entire curve could be drawn 
when 6 is known for any two values of N. It is preferable, however, to work 
out the largest possible number of values of 6 and draw an envelope to the 

plots of these values. After that the values of c and m could be determined 
by selecting two suitable points on that envelope. 

It is to be noted, however, that in hydrological phenomena the maximum 
positive deviations are likely to be relatively greater than negative deviations 
(see Art. 7). As far as the range (for a draft equal to the mean) is concerned, 
we should adopt the absolute maximum whether it is positive or negative. 

For the case of a draft less than the mean, on the other hand, negative devia- 
tions only are significant. For practical purposes, therefore, the deviation- 
curve should be constructed with reference to negative deviations only. 


Application of the Deviation-Curve 


In the ideal case of a draft exactly equal to the mean, the range may be 
determined graphically as shown in Fig. 7. Given the deviation-curve AB, 
we can draw another curve CD representing the product Né. By the construc- 
tion shown in the figure, a third curve, EF, representing the ratio R/Mg can 
be drawn, whence R is determined. It is understood that, in the absence of 
knowledge of the absolute mean, the overall mean for the whole of the existing 
‘record would have to be used instead. It is of interest to note that, despite 
this fact, the deviation-curve will automatically indicate the maximum devia- 
tion for that mean. 

If 6 is expressed as in (16), 


wid] = N-S =e 
n@ 


which is the equation of curve CD. The equation of curve EF will be: 


R ¢__ yi-@ (18) 


2 I-n 


a 
= 


For random events, m = 0.5 and, in consequence, R/M, varies as y N,as 


| 
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_ has already been shown. For natural phenomena, on the other hand, m would 
generally have some value between 0 and unity. The former value would 
signify that R varies linearly with N and the latter would Signify that R has a 
limiting value independent of number (as in the case of a perfectly periodic 
occurrence), 

Since 6 cannot possibly increase with N, m cannot be negative and, ac- 
cording to (18), R cannot vary as N raised to a higher power than unity. The 
fact that in Table III such powers are encountered is due to inconstancy of the 
standard deviation. It is clear that if o varied progressively with N, the ex- 
ponent of N in (4) would be affected. 

As regards the case of a draft less than the mean, it has already been 
pointed out that there is no justification for supposing the existence of a fixed 
relation between the storage needed in that case and the ideal range. It may 
be further noted that the crucial question in that case would not be the maxi- 
mum deficit within a given period but the maximum deficit with respect to 
the reduced draft B that might be encountered at any time. 


In general we have: 


S = Ng (B - Mg) (19) 


where Ng is the length of the maximum-deficit period and Mg is the corres- 
ponding mean. This relation is demonstrated schematically in Fig. 8. Put- 
ting B = relative departure of B from Mg and $g = relative departure of Mg, 


(19) reduces to: 


5 
Fyre (53- B) (20) 


Evidently, the value of S or S/Mg we should seek is that for which (N 6 - 
NB) is a maximum. Now, according to (17), 
WSs. ¢. wo 
Differentiating (20) with respect to N and equating to 0, we get: 
| e(l-m)n™ = B (21) 


whence the critical value of Ng would be: 


£ [ stam ] tm (22) 
Ng = B 


By inserting that value in (20), the required value of s will be obtained. 
The graphical solution of that problem is shown in Fig. 9. The (N6)-curve 
in that figure is the same as that given in Fig. te Ss straight line OG rep- 
‘resents the quantity (NB). The critical value of S/M, is given by the maxi- 
“mum intercept between the two plots. ss 

If m is very small or zero, Ng will be very large or infinite. In sucha 
“case, the choice of a suitable value of S will have to be based solely on prac- 


tical considerations. 


ry 


* 


1082-20 HY 5 October, 1956 


Ojo 200 Sor stor Som 60 747 ao. ae ce 
(N) 
Fig. 9 — DERIVATION OF THE M/XIMUM DEFICIT 


(2) 


Fig. I0 = REGULATION FOR A DRAFT LESS 
THAN THE MEAN 


ee ene eR a a) 


ASCE FATHY - SHUKRY 1082-21 
Program of Regulation 


Having now come to the conclusion that the draft should be fixed always at 
some value smaller than the probable future mean, the program of regulation 
will be as illustrated in Fig. 10. The full line Oabcd represents the mass- 
curve for the supply and the broken line Oa'b'cd' represents the mass-curve 
for the actual draft. During such periods as ab or de, the reservoir will be 
full and the whole of the natural supply will have to be passed downstream. 
Evidently, regulation cannot be expected to proceed according to plan until 
the reservoir has been once completely filled. A special program will, of 
course, have to be laid down for the initial filling. 


Case of Variable Draft 


In some cases, the draft from a reservoir may have to be varied from 
year to year instead of being maintained at a constant figure. A case in point 
is that of the proposed Lake Albert Reservoir on the Upper Nile. It is known 
that the Main Nile receives its waters from a number of important sources 
besides the Equatorial Lakes, namely, the Bahr el-Ghazal, the Sobat, the 

Blue Nile and the Atbara. With the exception of a small fraction of the Blue 
Nile supply (that coming from Lake Tana) the flow from these sources cannot 
be put under over-year control.* As far as Egypt’s requirements are con- 
cerned, there would be no sense in maintaining a constant outflow from Lake 
Albert. The function of the Albert Reservoir should be to make up any de- 
ficiency in the supply reaching Egypt from other sources below the require- 
ments. Thus, the outflow from that reservoir would have to be varied be- 
tween wide limits. 

To determine the necessary capacity in such a case, we need a sample ac- 
count of the annual requirements from the reservoir (compiled from past re- 
cords of the supply from other sources) in addition to the record of the na- 
tural inflow. We should see to it, by making a proportionate reduction in the 
demand if necessary, that a safe margin is left between the mean draft and 
the mean-supply. 

One way of dealing with that problem is to compute the departures in the 
draft from its mean, transfer these departures with reversed sign to the sup- 
ply and construct the deviation-curve for the latter on the basis of the modi- 
fied figures. The maximum deficit could then be determined as explained in 
Art. 12, with reference to the mean draft. 

Another way is to construct a deviation-curve for the draft separately, re- 
duce the scale of 6 in the ratio B/Mg, where B = mean draft, and plot the 
resultant (N6)-curve below the mean-draft line, as shown in Fig. 11. The 
value of S/M, required will be the greatest intercept between that curve and 
the (N6)-curve for the natural supply. If a maximum is not found within a 
period of, say, 50 years, then the intercept at N = 50 may be adopted. 

_ The second method is the more safe as it gives the worst possible combi- 
nation of an excess in the demand with a deficit in the supply. The probability 
of such a combination may be low when the supplies from the different 
sources are not correlated. The question, however, is one of comparison be- 
tween the additional capital outlay to be incurred and the damage that might 
ensue in the event of that coincidence actually taking place. 


* Unless the Aswan High-Dam Scheme, which is at present under considera- 
_ tion, proves feasible. In that case, the Albert Reservoir will not be needed. 


: 
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Effect of Sub- Annual Fluctuations 


As far as hydrological phenomena are concerned, the natural time unit of 
an individual observation in the consideration of long-term storage is one 
year. The theoretical treatment of the problem, however, would not be com- 
plete without considering the probable effects of sub-annual fluctuations on 
the result. It is not difficult to see that shortening of the time-unit might 
lead to an increase in the range of cumulated departures. 

In this respect, it is to be noted that while cyclic variation in annual rain- : 
fall or run-off observations is generally erratic, practically all such phenom- 
ena are subject to more or less definite sub-annual cycles, That calls for an 
independent treatment of sub-annual fluctuations in either the supply or the : 
demand. 

In a given reservoir, the effects of sub-annual fluctuations are likely to be 
felt only when the contents are about to touch either their minimum or their 
maximum value as indicated by the yearly account. Some additional capacity 
may be needed to allow for secondary deviations within the year preceding 
either event. That capacity can be determined on similar lines to those fol- 
lowed in the design of sub-annual reservoirs. 

In the majority of cases, however, the need for such a correction may be 
eliminated by a suitable choice of the date at which the hydraulic year is sup- 
posed to begin, If that date is so chosen that there would be no secondary 
filling during the last few months of the year, the reservoir contents would be 
most unlikely to fall below their ultimate value within the year at the end of 
which they are supposed to drop to their minimum. Whatever happens at the 
top end does not matter much since there would always be some means of 
disposing of excess water. 

It is to be borne in mind, however, that in any work designed on probability 
considerations, matters cannot be expected to run exactly according to 
schedule. Thus, in the running of an over-year reservoir, the draft should be 
modified according to some sliding scale as the contents approach their pre- 
vised minimum value. That might render a highly refined theoretical treat- 
ment unnecessary. 


Flood-Storage and Escapage 


Since, in practice, the mean draft adopted should always be less than the © 
predicted mean supply, the reservoir would, as a rule, have to remain for 
long periods on end at or near its maximum level. That is clearly demon- 
strated in Fig. 10. In a small reservoir, any sudden rise in the inflow during 
such periods might be met by instantaneous escapage through a spillway. In 
a reservoir of large surface area, the provision of a spillway of sufficient 
evacuating capacity to prevent an appreciable rise in the water level would — 
not be practical. It might also be inadvisable to install an automatic spillwa y 
if the sudden release of large volumes of water from the reservoir is likely — 
to cause damage in the downstream channel. # 

Thus, any large -capacity reservoir would, as a rule, have to be worked is 
asa flood escape in addition to its normal function. Knowing the limiting 
value of the discharge that could be safely passed downstream and the maxi- 
mum flood discharges (into the reservoir) to be expected, the computation 
the flood-storage needed is a simple matter. For obvious reasons, the pro- 
gram of flood regulation should be so arranged that the relief storave of a 
flood be disposed of completely before the arrival of the next flood. 
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(NV) 
Fig. II = STORAGE FOR A VARIABLE DRAPT 


(S) - 


(N) Years 


Fig. 12 = DEVIATION-CURVE FOR THE NILE AT ASWAN. 
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Effect of Storage Losses 


When a long-term reservoir is designed so as to guarantee a certain 
draft, that draft would obviously include storage losses. We may write: 


1- p= 224 (23) 


where B here denotes the useful draft and L is the rate of loss at the adopted 
storage level. Hence, 


B = ¥.(z= Pp) = 25 (24) 


It is clear that the closer the gross draft approaches to the mean supply, 
the greater will be the storage capacity needed. Since, however, storage 
losses also increase with the capacity, it will generally be found that there is 
an optimum value of the capacity beyond which B would begin to decrease 
again. That feature is clearly demonstrated in the practical example given 
in the next section. 


Engineering Application 


Factors Governing Design 


It has been said before that in the preparation of a long-term storage pro- 
ject, the principal tools in the hands of the designer are sound practical judg- 
ment and a thorough understanding of the principles involved. In the forego- 
ing treatment of the problem, the relevant physical factors have been 
discussed at some length in order to make these principles as clear as pos- 
sible and to reveal the considerations in respect of which the designer has to 
use his judgment. At first glance, the method of treatment suggested might 
appear to be time consuming. Actually, it will be found to be more rapid than 
any un-guided trial and error method, in addition to giving the designer a 
surer footing on which to stand. 

The keynote of the treatment is the Deviation-Curve. That curve gives us 
as complete an account as possible of the physical characteristics of the phe- 
nomenon under consideration. It is admitted that great care is needed.in the 
preparation of that curve since it is the basic factor in the whole problem. _ 

Once the deviation-curve has been prepared, the capacity-yield relation- © 
ship for a wide range of both can be readily obtained. From relations (18), 
(19) and (20) of Art. 12 we get: 


1-m 
ag oi tari or 4 
My em | 6 | (25) : 


By assigning different values to 8, the corresponding values of S/Ma, ors 
(since Mg is known), can be readily computed. 

The relation between S and the rate of storage loss L must, of course, be. 
known. A mathematical expression is not necessary, since the problem can — 
henceforth be solved graphically. With the aid of (24), we can draw a curve — 
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representing the variation of the effective draft B with S. On the basis of 
that curve, another one representing the cost-yield relationship can be ob- 
tained. The choice of capacity may then be made in full view of all the physi- 
cal and economic factors involved. 


Practical Example 


As the present investigation was originally inspired by schemes for over- 
year storage in the Nile Basin, it would be natural to use one of these 
schemes as an illustration. The famous High-Dam scheme is well suited for 
our purpose. The working storage arrived at through use of the deviation- 
curve method, as will be seen presently, is sensibly smaller than that 
adopted in the officially approved design of that scheme. 

The scheme involves the creation of an over-year reservoir in the valley 
of the Nile just upstream of the present Aswan Dam. The whole of the flow 
of the River, including the silt-laden flood water, will be impounded in that 
reservoir. To allow for loss of capacity through silt deposition, a “dead” 
storage of 30 milliard cubic metres is to be provided (one milliard = 109). 
The “live” or working storage, according to official accounts, is to be 70 
milliards m$. On top of that, there will be added another 30 milliards for 
flood relief, bringing up the gross capacity to 130 milliard m3. 

_ The mean natural-river supply at Aswan for the period 1871-1950 is 93 
milld.m3/ annum. Deducting 4 milliards against Sudan abstractions, the most 
probable value of the mean supply during the next hundred years or so will 
5e 89 milld.m3/annum. If conservation works are carried out in the Upper 
Nile Basin, as projected, these conditions will change. That possibility may 
nave been allowed for in the actual design. In the present example it will be 
ignored, and the figure of 89 millds. will be supposed to represent the abso- 
lute mean supply (Mg). 

Estimates of reservoir areas and capacities at different levels have been 
orepared by the High-Dam Authority. Figures have also been given for stor- 
age losses at different levels. Evaporation losses are not difficult to esti- 
mate but it is impossible, at present, to make accurate forecasts of seepage 
ind saturation losses. Evaporation in the Aswan district is known to be about 
1000 mms. per annum. In the following computation it will be assumed that 
he gross rate of loss will be equivalent to a depth of 3 metres per annum on 
he water surface. Precipitation in the district is practically nil. 

In the computation of the storage losses for different reservoir capacities, 
wo special corrections have had to be made. Firstly, since at any adopted 
‘apacity the reservoir would not remain always at its highest level, some re- 
juction in the corresponding rate of loss should be allowed for. Secondly, 
since the working storage is to be superposed on the dead storage referred to 
ibove, the corresponding loss should be reckoned on a figure 30 milld.m3 
rreater. 

_ The solution of the problem in accordance with the foregoing data is given 
‘raphically in Figs 12 and 13. Fig. 12 shows the deviation-curve for the 
atural supply at Aswan after correction for Sudan abstractions as mentioned 
bove. The observed maximum deviations are well covered by the adopted 
urve of which the equation is: 


- rw (26) 
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0-85 


Pig. 13 - YIELD-CAPACITY CHART FOR THE ASWAN HBIGH-DAM, 
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‘Thus, according to the notation used, c = 0.4 and m = 0.333. 
In Fig. 13, curve (ab) is a plot of the quantity (1-8), which is equivalent to 

(B+L)/Ma, against S/Mg. Curve (cd) is a plot of L/Mg on the same base. 
Curve (ef) represents the effective-draft ratio B/Mg. In the last plot we see 
vat once that between the values of S/Mg: 0.5 and 1.0, the gain in effective 
draft decreases rapidly and, thereafter, B remains practically constant des- 
pite the increase in S. It is clear that if the plots were extended B would 
€ventually decrease instead of increasing with S. The disposition and char- 
acter of the B/Mg-curve evidently depend very largely on the nature of the 
L/Mg-curve. That calls for very careful study of the question of storage 
losses. 

No data regarding the cost increments for different storages are at hand 
but it would appear from the figure that an increase in the value of S/Mag 
above about 0.5 would not be economical, At that value, the useful draft is 
0.791 Ma, or 70 milliards m3/annum and the working storage is 45 milliards. 
By raising the storage to 70 milliards (0.785 M,) the gain in draft would be 
Only about one milliard per annum. 

The final result may be summed up in the following: 


a) S = 45 milliard m3 (gross storage = 105 millds.) 
b) B= 70 milld.m3/annum. 
c) B = 0.138 
d) Ng = 7.2 years. 
The last figure shows that deviations for short periods may be much more 
important than deviations for long periods. In the construction of the 
deviation-curve, therefore, particular attention should be paid to short-period 


deviations. 

In conclusion, it may be pointed out that the sensitivity of the result to the 
physical characteristics of the phenomenon considered, as may be gathered 
from the above example, renders the use of empirical capacity-formulae of 
a general character inadvisable. 
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Discussion of 
“TIDAL COMPUTATIONS IN SHALLOW WATER” 


by J. J. Dronkers and J. C. Schénfeld 
| (Proc, Paper 714) 


| J. J. DRONKERS! and J. C. SCHONFELD.2—In his comment, Mr. Fenwick 
gives an interesting supplementary note on the question of the true simulation 
of the flow pattern in a distorted model. His remark applies to relatively 
large-scale models with moderate distortion. 

In the writers’ comparative discussion on models and computations, they 
had more particularly in mind strongly distorted, small-scale models. The 
following consideration will make this clear: 

It is a useful principle in tidal hydraulics, to distinguish between the gener- 
al tidal motion as characterized by the fluctuating levels and total flows, and 
‘more local phenomena defining the detailed distribution of the total flow. 
These two aspects of the tidal motion need not necessarily be treated by the 
same means. 

In case of a tidal system of small extent, a hydraulic model on a large 
scale and without or with moderate distortion, may meet all requirements. 

_In case of complicated and extensive tidal systems, however, the demands 
_imposed upon the model from the two points of view, may be so divergent, 
that it is indicated to make a small-scale, strongly distorted model of the en- 
tire system for the investigation of the general tidal motion, together with one 
or more large-scale models of particular parts of the system for the investi- 
gation of the local problems. 

The great value of the larger-scale hydraulic models for the investigation 
of the local tidal problems, is beyond discussion. However, for the investiga- 
tion of the general tidal motion, the hydraulic model is not the only possible 
tool, and an electric analogue or a mathematical “model” may serve as a 

substitute for the small-scale hydraulic model, or as a supplementary check. 

As has been stated, the comparative discussion in the paper referred in 
particular to hydraulic models of large tidal systems, on a small horizontal 


scale. 


4 Chief Mathematician, Central Research Div., Netherlands Rijkswaterstaat. 
2. Chf. Engr., Central Research Div., Netherlands Rijkswaterstaat. 
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Discussion of 
“THE IMPORTANCE OF FLUVIAL MORPHOLOGY IN HYDRAULIC 
ENGINEERING” 


by E. W. Lane 
(Proc. Paper 745) 


E. W. LANE,! M. ASCE.—The author very much appreciates the discus- 
sion of his paper by Messrs. Kuiper, Nimmo and Happ, and believes that the 
points they have brought out have added materially to its value. 

Mr. Kuiper has emphasized the difficulties of making accurate estimates 
of the magnitude of the changes in streams of the nature discussed in this pa- 
per. It is true that in the present state of our knowledge, only under the most 
favorable conditions is it possible to make reasonably accurate quantitative 
estimates. However, by an application of the general principles discussed in 
this paper, it is frequently possible to determine positively that an undesir- 
able action will take place. When this is the case the engineer cannot ignore 
it. He must make the best estimate of the magnitude and rate of change that 
the present state of the science and data available will permit. Although the 
estimates in some cases may be far from accurate, the results on the whole 
will be much better than if the action was entirely ignored. The engineer 
must, of course, be perfectly frank regarding the degree of uncertainty of his 
estimates. Often very rough estimates can be very valuable. For example, 
there is the case of a very expensive recreational development which was 
built at the upper end of a reservoir lake, which had only mud flats in front of 
it a few years after it was built. There is another case of a similar situation 
in which great loss would have resulted had not the plans been abandoned. In 
either case the result could have been foretold with very limited quantitative 
studies, and it is not unlikely that in the second case the plan was abandoned 
because such an estimate was made. As recorded observation on actions of 
this type increases, and the science of sediment transportation advances, 
more and more accurate estimates will become possible. 

The writer feels that Mr. Kuiper’s example presents a more pessimistic 
Picture than occurs in many cases. The rate of aggradation fifty miles up- 
stream from a reservoir as posed in his problem, would depend largely on 
the rate at which the reservoir filled with sediment, and this would probably 
be controlled largely by the suspended sediment, the load of which often can 
be measured with reasonable accuracy. 

In other problems of the nature under consideration the bed material load 
is often the main factor in causing the changes. Its magnitude can frequently 
be estimated with reasonable accuracy from suspended sediment measure- 
ments, and a “rating curve” of this load for the various discharges obtained. 
By combining this curve with a discharge duration curve of the stream at that 
point, the average annual bed material load can be computed. Approximate 
Besics of the volume of deposition after a certain period can often be made 


i. Cons. Hydr. Engr., Fort Collins, Colo, 
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by drawing profiles representing the nature of the changes from the original 
profile which one would expect from the figures presented in the author’s 
paper. Computed backwater curves can be used to indicate the upper limit of 
the effect. The volumes enclosed between this profile and the original one 
can be computed, and the probable time required to make this change com- 
puted from the relation of this volume to the average annual load. The ac- 
curacy of such estimates depends to a large extent on the knowledge and ex- 
perience of the one making them, but they are likely to produce much more 
satisfactory results than would result from ignoring the changes. In drawing 
such an altered profile, a study of the changes which occurred in other cases 
is very helpful, and the author would like to emphasize Mr. Kuiper’s state- 
ment of the value of systematically collecting data on river regime and the 
effect of artificial interference with it as a means of perfection of the ac- 
curacy of future estimates of the magnitude of such effects. 

Mr. Kuiper has correctly pointed out the difficulty, under present condi- 
tions, of obtaining accurate results by the direct application of sediment 
transportation formulaeto these problems. The author feels, however, that 
if Mr. Kuiper presented in detail the results of the excellent work that was 
done in analyzing a problem of this general nature by Mr. Kuiper and col- 
leagues in connection with a reclamation project on the Saskatchewan River 
in Canada, the readers would agree that more reliable solutions were pos- 
sible in some cases than were indicated by Mr. Kuiper’s discussion. 

Mr. Nimmo has correctly called attention to the effect of attrition and 
selective deposition on the shape of stream profiles. The author did not in- 
clude this action in his paper since it would have unnecessarily complicated 
the presentation of the ideas he wished to introduce. The magnitude of this 
effect can usually be obtained from the shape of the stream profile in its in- 
itial condition, and the changes of the nature considered in this paper can be 
treated as changes from this original profile. 

Attrition and selective deposition are important aspects of geology. By 
observing the gradual decrease in particle size with distance in sediments 
laid down many millions of years ago, geologists can show that mountain 
ranges formerly existed where little evidence of them is present today. Ex- 
periments show that attrition is appreciable in the transport of sediment of 
large sizes, but is very small for material of sand size.* In a given stream 
it would be difficult to separate the effects of attrition and selective deposi- 
tion. In general it is probable that where the valley is narrow, so that the 
opportunity for storage of sediment is limited, the attrition is the predomi- — 
nant action, but where large storage spaces are available, selective deposi- — 
tion is the major effect. Attrition has been extensively studied** and : 


* C.C, Inglis and D. V. Joglekar—Rate of Sand Attrition in Channels, Indian 
Waterways Experiment Station, Poona Publication 8 pp 62-64, 1945, 

** S, Shulits—Fluvial Morphology in Terms of Slope, Abrasion and Bed Load- 
American Geophysical Union Transactions, Vol 17, Part 2, pp 440-444, 
July 1936. « 
S. Shulits—Rational Equation of River Bed Profile—American Geophysical 
Union Transactions, Vol. 22, Part 3, pp 622-630, 1941. 

WwW. C. Krumbein—The Effect of Abrasion on the Size, Shape and Roundness 
of Rock Fragments—Jour. Geol. Vol. 49 No. 5, pp 482- 520, July/Aug. 194 

A. O. Woodford—Stream Gradients and Monterey Sea Valley, Bull. Geol. ‘ 
Soc. of America, Vol. 62, pp 799-852, 1951. 
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Separation of the two effects in any case could probably best be made by 
computing the attrition effects and assuming that the remainder of the change 
was due to selective deposition. 

The picture and description of Coopers Creek in South West Queensland 
presented by Mr. Nimmo is very striking. It probably represents an extreme 
case of braiding in a channel, due to aggradation, where the spreading of the 
channels over a very large width was permitted, 

Dr. Happ’s discussion is valuable in illustrating the difference in the 
thinking of the geologist and the engineer due to the relatively minute time 
involved in the cases handled by the engineer when compared with those usu- 
ally considered by the geologist. He is right in calling attention to the fact 
that the classification of streams using the terms youth, mature and old age 
is only descriptive, but the use of these terms descriptively by engineers 
would be useful in many cases. Dr. Happ is also correct in noting that 
stream flow variation is a factor in many morphological problems encoun- 
tered by the engineer. The effect of variable flows as compared with uniform 
ones, is that the variable flows transport higher sediment loads for the same 
total quantity of water involved, than uniform flows. The effect of increasing 
or decreasing the variability of flow, would have the same effects in the 
stream action as increasing or decreasing respectively, the discharge of the 
stream, 

The author’s classification of the morphological changes of a stream was 
drawn up to present the cases which are likely to be encountered by the engi- 
neer in a way which would be readily grasped by him. Other classifications 
are possible, and for a pure science or geological approach, perhaps Dr. 
Happs classification would be better. His classification could be reduced to 
six by combining the two items under the same letter. For example, the two 
items designated B could be combined into one stated: change of slope caused 
by lowering or rise of the base level. The author’s six cases could similarly 
be reduced to three, 

_ For engineers, the principal objection to Dr. Happ’s classification is his 
inclusion diastrophic uplift or crustal tilting. The hydraulic engineer con- 
tinually encounters cases where crustal tilting which has taken place in the 
past is important, but the author knows of no case where the crustal tilting 
which may occur during the life of an engineering project has to be taken in- 
to account in its design. 

_ Dr. Happ has called attention to the use by geologists of a second concept 
of grading as applied to streams which differs fundamentally from the equili- 
brium concept discussed by the author. This second concept involves smooth- 
ing out of the stream profile and is not concerned with equilibrium or balance. 
The author’s treatment largely followed the balance idea, which is presented 
by Kesseli(21)* and Mackin,(25) The smoothing out concept of grading is 
probably not widely used since the author did not find it in his rather exten- 
sive search of geological literature. Since reading Dr. Happ’s discussion he 
has encountered a possible case of it in the definition of “graded” by L. 

La Forge given in A Glossary of the Mining Mineral Industry.** 

7 Graded is defined by La Forge as “Brought to or established at grade, 
through the action of running water carrying a load of sediment, by eroding 
br degrading at some places and depositing or agegrading in other places.” 


2 These references refer to the original list. 


r U. S. Bureau of Mines Bulletin No. 95—A. H. Fay. 
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This eroding at some places and depositing at others may have involved in 
La Forge’s mind a smoothing out of the profile, but this is questionable since 
his definition of “grade” is “That slope of the bed of a stream, or surface 
over which water flows, upon which the current can transport its load without 
eroding or depositing.” 

In view of the two different concepts of grading held by geologists, Dr. 
Happ is probably correct in suggesting that engineers use the term “poised” 
suggested by G. H. Mattlies Hon. MASCE(16) rather than graded, when desig- 
nating a stream in balance. Another term to describe this condition which 
may be useful is “quasi-equilibrium,” as used by Leopold and Miller.* How- 
ever, regardless of what is the best terminology, the engineer is indebted to 
geologists for much valuable work in this field. : 

Since writing the original paper, it has occurred to the author that a better 
way than the “equation” Qod ~~ QyS to present the relations between the fac- | 
tors involved would be to use an analogy to a balance, in one pan of which was’ 
Qs and d and in the other pan Qw and S, as shown in Fig. 1. To reestablish an 
equilibrium which existed in a stream between these four factors Qs, d, Qw 
and S when it was upset by a change in any one of them, it would be necessary 
to make a change in the opposite direction inthe factor in the same pan as the 
changed factor causing the disequilibrium or change in the same direction as 
the disturbing factor in one or both of the factors in the other pan, or a com- 
bination of these three changes. These changes are of the same pattern as 
would have to be made in the analogous case of a balance in Fig. 1 if the vari- 
ables were considered to represent weights. For example, an equilibrium 
which was upset by an increase in Qs could be restored by a decrease in the 
values of the factor d which is in the same balance pan or an increase in 
either or both the factors Qw or S which are in the other pan, or a combina- 
tion of decrease in d and increase in Qw and/or S. 


Fig. 1 
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* L.B. Leopold and J. P. Miller—Ephemeral Streams—Hydraulic Factors 
and their Relation to the Drainage Net, U. S. Geological Survey Profes- 
sional Paper 282-A, 1956. 
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Discussion of 
“MINIMUM PRESSURES IN RECTANGULAR BENDS” 


by M. B. McPherson and H. S. Strausser 
(Proc. Paper 747) 


M. B. MC PHERSON,! A.M. ASCE and H. S. STRAUSSER,2 J.M. ASCE.— 
The authors are indebted to Professors Shukry, Robertson and Rice for the 
supplemental information and clarification given in their discussions. Their 
unanimous endorsement of the method presented for predicting minimum 
pressures was particularly heartening. 

It had been hoped that the discussion would be limited to closed rectangu- 
lar bends. However, Professor Shukry presented some valuable information 
related to subcritical open channel flow, and Professor Robertson presented 
hitherto unavailable limiting Reynolds number data for elbow meters. The 
Suggestion by Professor Rice on the use of vanes to minimize cavitation de- 
serves serious consideration by both designers and researchers alike. 

Concerning effects of viscosity, the lowest Reynolds number in tests per- 
formed on the Waynesboro and Mt. Alto models was approximately 105, Be- 
cause of limitations in size or velocity dictated by components other than 
elosed bends in model studies of hydraulic structures, it appears that data for 
values less than 10° would seldom be encountered. If the data of Figure A by 
Robertson is applicable as well to closed bends of rectangular cross-section, 
the lowest value for which such model tests on bends should be run would ob- 
viously be 109, Unfortunately, the authors do not have substantiating ee 
The data of Figure 2 by Shukry are all for Reynolds numbers less than 10°, 

The paper was directed towards design rather than analysis. However, 
the corrections detailed by Shukry might provide a better appraisal of model 
results, particularly for low values of x. 

Robertson has emphasized the point that correlation with irrotational flow 
piezometric head, as presented, is limited to the station of the maximum 
piezometric head differential. Put another way, it would appear that the 
general curvature of the streamlines in a closed rectangular bend are closest 
to the curvature of the bend itself at a station about 45° from the P.C. Some 
streamline curvature of necessity extends beyond the bend curvature, and at 
every station other than that of closest coincidence, is less than the bend 
curvature. This apparent coincidence of streamline curvature at the station 
of lowest head is the only justification for using a frictionless approach in 
evaluating minimum pressures. 

Irrotational streamline curvature, as reflected in corresponding differ- 
ences in head, may be obtained analytically, graphically or by analogy. The 
over-all ae: analytical solution is beyond the mathematical capacity of most 
engineers, including the authors. Graphical solution by means of a two- 
dimensional flow net is entirely too inexact, particularly for bends of short 
radii. At the conclusion of some studies on irrotational flow through bends of 
€ 
1. Associate Prof. of Civ. Eng., Lehigh Univ., Bethlehem, Pa. 
. Asst. Prof. of Civ. Eng., Univ. of Washington, Seattle, Wash. 
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circular cross-section by the method of electrolytic-tank analogy, Murthy (3) 
recently performed a similar test on the Waynesboro model bend. The plas- 
tic bend, with one side removed, was filled to mid-width with electrolyte, and 
electrical potential drops were measured, The data obtained for the bounda- 
ries appears in Figure 5a, while the apparatus used is shown in Figure 6, 
since it may be of interest to some readers. For the analogy data, hy is the 
equivalent piezometric head before or after the bend in the zones of uniform 
(parallel streamline) flow. The rearranged data from Figure 5 of the paper 
is also shown, for comparison, Inasmuch as the total energy in the latter 
case decreases with distance, an equitable comparison is obtained by using a 
varying reference; the term he is the piezometric head as previously recom- 
mended for use; that is, V2/2g below the total energy line shown dashed in 
Figure 5. Note that any irregularities present in the hydraulic model were 
reproduced per se in the analogue. The value of Cx, obtained from the ana- 
logue was 1.50 vs. 1.45 for Figure 5 (1.415 in Table 1 of the paper was the 
average value of several runs with a variety of approach conditions) and 1.43 
from Equation 2. Despite probable inaccuracies due to construction and mea- 
surement, a reasonable comparison is obtained. Considering the inherent 
secondary motion, varying boundary shear and probable separation, it is in- 
deed surprising that there is any semblance of correlation to irrotational 
motion. According to Figure 5a, the piezometric head distribution for a 
closed rectangular bend is quite similar to that for irrotational flow. The ar- 
gument offered by Robertson concerning the proper limited use of potential 
theory is valid only with regard to a station having similar streamline and 
boundary curvature. The “free-vortex” is a limiting case with streamlines 
forming concentric circles. 

In design, the use of Equation A set forth by Robertson yields lower values 
of Cy, and hence lesser anticipated minimum pressures than Equation 2 of the 
paper, regardless of simplicity. Hence Equation A cannot be regarded as 
conservative for use with rectangular bends, especially for low values of x. 

Since 90° circular bends have been introduced to the discussion, some at- 
tention must be given to the differences in results obtained for both types of 
bends. The available data presented in Figure 2 of the paper shows both con- 
sistency and correlation with Equation 2 for the 45° station of rectangular . 
bends of 90° and 180° central angle. Similar data restricted to 90° circular — 
bends at the 45° station (Figure 2 of Reference 2 of the paper) show consider- 
able scatter, with representative data as often approaching Equation 2 as A, _ 
and with some data yielding values of Cy, well in excess of those for Equation 
2. From this limited evidence, the use of Equation 2 for predicting safe cir- 
cular bend performance could hardly be considered conservative, much less ; 
Equation A. The following Table includes data from two elbow meters tested 
with water at Lehigh, together with data by Murthy from half-section electri- 
cal analogy measurements on the same elbows (these two elbows were used 
since flow data gave unusually high values of Cy for given x than indicated by 
most of available data): 


3. “Potential Flow in 90° Bends by Electrical Analogy,” by D. S. N. pct: 
May, 1956, part of M.S. program, Lehigh University Library. 
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Tabulation of Values of Ck at 450° 


for 90° Circular Elbows 


x = 3.50 x - 4,50 
From flow test - 1.58 1.00 
By Murthy, 
electrical analogy - 1.55 1.02 
Using Equation 2 - 1.28 0.95 
Using Equation A - 1.14 0.89 


The data obtained by electrical analogy, above, should be a reasonable ap- 
proximation of the irrotational case. Only three bends (the highest values of 
x) out of six tested at Lehigh gave values surpassing those for Equation 2. 
The Cy values for all six exceeded Equation A, with x = 2.25 to 4.5. There- 
‘ore, with circular bends, the use of Equation 2 for design calculations does 
aot appear to be conservative. Further, it is evident that a corresponding 
2quation based upon irrotational flow in a circular bend would differ appreci- 
ably from Equation 2 and even more from Equation A. The mathematical de- 
-ermination of the irrotational field in a bend of circular section is more dif- 
‘icult than that for a rectangular section, if it is solvable at all. As far as 
results are concerned, it would appear that the points of separation in a 
rectangular bend must be relatively fixed, or close te the 45° station, The 
variety of results obtained with circular bends, on the other hand, might con- 
seivably arise from a more variable point of separation, or some other factor 
10t readily discernible. The authors have observed a relatively greater ef- 
‘ect on the head distribution as the result of upstream disturbances on circu- 
ar bends than on rectangular bends. 

From the available evidence, the design criteria appears to be more defi- 
1ite for rectangular bends, due to an apparent lack of sensitivity. 

The experiments by Mr. Murthy were sponsored by the Civil Engineering 
Department and conducted in Fritz Engineering Laboratory, Lehigh 
University. 
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RECTANGULAR BEND - WAYNESBORO 
6=90° 
e=Flow Test D-I-D 


Curve = Electrical Analogy oe eee Te 
Data by Murthy 
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ELECTRICAL ANALOGY APPARATUS 
(The Waynesboro bend is in the foreground) 


Figure 6 
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Discussion of 
“RAINFALL DEPTH-DURATION RELATIONSHIPS” 


by Herbert M. Corn 
(Proc. Paper 840) 


HERBERT M. CORN,! A.M. ASCE.—The discussions by Messrs. Hersh- 
field and Wilson primarily question whether the storms of the United States 
(Yarnell’s tabulation of the most intense storms),2 which were used in the 
‘Studies of rainfall regimes of French Morocco, were analogous and repre- 
Sentative. Further, questions are raised regarding the basis of extending the 
‘mean percental depth-duration curve, derived from the studies, beyond the 
Tange of data used for the formulation. Actually, both questions are perti- 
nent to possible future practical utilization of the curve, but neither ques- 
tionable practice was applied directly to the problem nor intended for use in 
that manner, without proper modification for regional differences. It is 
fundamental that areas with a high incidence of thunderstorms or having low 
annual precipitation totals, whichever is characteristic, would not be used to 
determine regimes of other meteorologically dissimilar areas. 

As noted in the subject paper, the studies described therein were con- 
cerned principally with determining the existence of probabilities in the logi- 
cal appearance and sequence of incremental maximum subintensities within 
rainfalls, regardless of physical or geographical origin. It was desirous to 
apply the conclusions derived from the studies to the problem at hand, when 
appropriate methods for their application had been established. This, in fact, 
was accomplished, in Morocco, using rainfall records collected by the local 
weather observation stations and government experimental farms. Unfortu- 
nately, the background data, used at that time for the analysis, is not in the 
possession of the author for presentation at the present, however, the results 
therefrom do compare favorably with later frequency-intensity-duration an- 
alysis prepared by the Cooperative Studies Section, Hydrologic Services 
Division, Weather Bureau, U. S. Department of Commerce.? The latter 
studies indicated the 2-year one-hour rainfall, for two particular airbase 
locations in Morocco, to be 0.60 inches for one and 0.50 inches fox the other. 
The 2-year one-hour rainfalls that were determined for the same locations, 
by utilizing 6-hour local data and interpolating therefrom by use of the curve 
developed during the studies of short-time precipitation from stations in the 
United States, were 0.80 inches and 0.65 inches respectively. The depths 
deduced by both studies are favorably comparable and certainly within the 
realm of practical accuracy. 
| 


. Project Engr., Porter-Urquhart, McCreary & O’Brien, Cons. Engrs., 

_ Newark, N. J. at 

“Rainfall Intensity-Frequency Data,” by D. L. Yarnell, Misc. Publication 

_ No. 204, U.S.D.A., 1935, pp. 9-23. 

. “Rainfall Intensities For Local Drainage Design in Coastal Regions of 
North Africa, Longitude 11° W to 14° E. For Durations of 5 to 240 minutes 
and 2-, 5-, and 10-year Return Periods,” U. S. Weather Bureau, Dept. of 


Commerce, Washington, D. C., 1954. 
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The author is engaged, at the present time, in hydrological studies in con- 
junction with a highway design program in Honduras, Central America. In- 
tensity-duration-frequency curves are being prepared for the project by 
utilizing the methods of synthetic development previously used successfully 
in North Africa. It is proposed to describe the methodology and submit the 
background data and analysis to ASCE for publication upon its completion in 
the near future. 

The assumption, by the discussors, that the storms selected from Mr. 
Yarnell’s studies covered only the region of Central United States is not en- 
tirely correct. Actually, the data ‘used, for the analysis, were the first sixty 
most-intense rain storms tabulated in the publication2 and they are repre- 
sentative of geographical areas covering thirty-six states from Eureka, 
California to Eastport, Maine and from Baker, Oregon to Apalachicola, 
Florida. Therefore, the thunderstorm-day, dewpoint, or 2-year, 24-hour 
precipitation comparisons of Central United States and Morocco as noted in 
the discussions are inapplicable. Furthermore, as previously mentioned, it 
was realized that modifications of the conclusions derived from the studies 
would be necessary to compensate for the differences in climatic conditions, 
when methods for their practical application had been determined. 

The sixty storms analyzed were, in fact, generally convective and mostly 
summer rains. Although, the selection of this type of rainfall was coinciden- 
tal with its position in the tabulation by Mr. Yarnell, it was of particular in- 
terest with regard to the engineering design of drainage facilities in Morocco. 

The atmospheric circulation affecting the climate of North Africa is pri- 
marily the west-east drift of the Atlantic Polar Front, which is generally lo- 
cated along the fortieth parallel in the winter. During this season, for most 
parts of Morocco, the greatest portion of the annual precipitation is accumu- 
lated. According to available meteorological information,4 these frontal 
rains are rarely of long duration, being on the fringes of passing cold fronts, 
and usually fall in the form of showers. Torrential rains are generally not 
observed during the winter and the incidence of thunderstorms, in that season, 
is rare. However, during the summer and autumn months thunderstorms are 
the most prevalent form of precipitation, in the areas studied, and average 
generally fifteen to thirty a year. The area, wherein the incidence of thunder- 
storms approached the higher figure, also, was characterized by a minimum ~ 
number of annual precipitative days; equalling on the average less than fifty. 
It follows, therefore, that the heavy rains resulting from thunderstorm ac- 
tivity, which often produce maximum runoff, would be of special interest to _ 
the engineer. Also, for this reason, it may be stated that a relevancy existed 
- between the studies of the convective storms from Mr. Yarnell’s tabulations — 
and the critical storms that could be expected to occur over the project areas, 

Another question raised by Messrs. Hershfield and Wilson, as it is under- 
stood by the author, was the advisability and validity of extending the mean r. 
percental curve derived from short-time precipitation data to cover longer 
periods, such as 12- or 24-hours. The author agrees that there is little basis 
for the extension of the curve beyond the limits of the foundation data, al- @& 
though in some cases where the application had been made, the resultant fit 3 
was, for all practical purposes, satisfactory, as shown in Figure 6 of the pa-_ 
per. This may be explained by the “within burst” distribution of rainfall 


4, “Global Physics and Meteorology in Morocco, Report on Our Knowledge i 
1947,” by G. Bidault and J. Debrach, 1947. 
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increments, as noted by the discussors. Similar relationships were found in 
an analysis? of long duration continuous rainfalls, although obviously the re- 
lationships will not apply to “among burst” distribution. At the present time, 
the studies of more than 160 long-duration storms, for time periods of from 
6-hours to 6-days, have been completed which conclude that similar proba- 
bility relationships exist for the rainfalls examined although different quanti- 
tative curve definitions were obtained, as had been expected. It is believed 
these studies corroborate the conclusions of the studies of short-time pre- 
cipitation described in the paper. 

No assumption was made by the author regarding the representation of the 
conclusions from the analysis of short-time precipitation to 24-hour, or 
longer rainfall regimes of remote areas. The discussors’ comments on the 
periods of “no-rainfall” for a normal 24-hour time period is correct, but to 
carry the thought to an extreme, comparisons of this type could be made be- 
tween one-day, monthly and annual rainfalls, which would be unrealistic. 

The author is aware of previous studies made by Messrs. Jarvis® and Lins- 
ley? and the Corps of Engineers regarding the uniformity of percental rela- 
tionships of storm rainfall and annual totals for certain areas, but did not in- 
tend to imply the generalization of the mean curve derived from the studies 
of short-time precipitation to an application of this magnitude. 

The intent of the curves of Figures 3 and 4 of the paper have been ap- 
parently misunderstood, somewhat, by the discussors with reference to the 
similarity and the 1- to 24-hour rainfalls, although the ratio has general ap- 
plication. The curves of Figures 3 and 4 indicate that for continuous (“within 
purst”) precipitation the ratio of the maximum 5-minutes to a total 120- 
minute storm will be the same as the maximum 1-hour to a total 24-hour 
storm. This general application is evident from examination of Table 2 and 
Figure 6 of the paper. The comparisons of 1- to 24-hour rainfalls of Central 
and the entire United States, which were cited, is predicated upon regional 
tlimatic differences of the areas and the lack of continuity of precipitation, 
which is normal to the longer periods as noted in the discussion. 

Figure 8, of the discussion, showed the regional differences in the 10- and 
30-minute rainfalls, between Zones I and II, of western United States. The 
mean ratio of the 10-minute rainfalls to 50-minute rainfalls, based upon the 
summating and averaging of ratios for both zones, is approximately 0.36. 

The ratio for these periods as determined by the curve shown in Figure 4, of 
subject paper, is 0.35. Further, the mean ratio for the period in question, 
1s shown by the curves of Figure 9 of the discussion which had been men- 
ioned in the paper for reference® is approximately 0.38. Consequently, it is 
elieved the comparisons herein noted are reasonable evidence of the proba- 
vility in the uniformity of intensity-duration relationships expounded in the 


aper. 


. “A Study of Rainfall Depth-Duration Relationships,” unpublished thesis for 
_ Master of Engineering Degree by Capt. C. J. Cox, Agric. and Mech. Col- 
lege of Texas, May, 1955. 

. “Rainfall Characteristics and Their Relation to Soils and Runoff,” Trans. 


_A.S.C.E., Vol. 95, 1931, pp. 379-423. 


,; “Frequency and Seasonal Distribution of Precipitation over Large Areas,” 


Trans. A.G.U., Vol. 28, No. 3, June 1947, pp. 445-450. 
. G. R. Williams, Hydrology, Engineering Hydraulics, P. 275. 
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The illustration by the discussors of the limitations of extrapolation of the 
curve beyond the range of data is noteworthy and is an excellent example of a 
common error committed by the uninitiated. The misuse of empirical formu- 
lae is commonplace and the author concurs that conclusions derived in this 
manner may not only be misleading but often can result in extravagance of 
design or worse, the failure of hydraulic structures. 

Unfortunately the regions, where the synthesis of maximum incremental 
periods would be appropriate, are typified by lack of detailed rainfall data. 
However, it has been the experience of the author that the technicians main- 
taining the observation stations in these areas are rarely unable, after years 
of experience recording data, to satisfactorily bracket the durations of par- 
ticular rainfalls by mannerisms of notations or remarks noted in the logs. 
For instance, in Honduras the accumulated rainfall depths are logged at the 
end of each hour and in most cases notations are made regarding the esti- 
mated portion of the preceding hour within which the rain occurred. This is 
particularly true of heavier rainfalls. Obviously, information of this type 
cannot approach the accuracy of mechanically recorded data. However, until 
the latter becomes available, a means of deriving intensity-duration relation- 
ships is possible and practical for use in these areas. 

The dissimilarity of measurements for stations of rather close proximity 
and the inconsistencies of data for the same station was a condition that 
generally existed in most locations from which data was available and ex- 
amined in Morocco. The constancy of its appearance precluded the possi- 
bilities of their being associated with convective precipitation only. It is 
noteworthy to mention that the Cooperative Studies Section, Hydrologic Ser- 
vices Division, Weather Bureau, U. S. Department of Commerce also found 
this condition to be troublesome in their studies in the same region. 

The storm of June 21, 1933, at Burlington, Vermont, was considered to be 
a peculiarity, by the author, inasmuch as the behaviour of the maximum in- 
cremental periods, of this particular storm, was entirely different than the © 
patterns of all the other storms used in the analysis and comparisons there- 
after. The fact that storms of these types regularly occur at Burlington does 
not increase its weighted value in a probability analysis when compared to 
the weighted averages and normal similarities of all the other storms ex- 
amined, 5, 

The author is grateful and appreciative for the response and comments of 
the discussors and hopes the paper and discussions serve as a small step in 
the direction of closer association between the meteorologist and hydrologist. 
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| Discussion of 
: “RESEARCH NEEDS IN SEDIMENT HYDRAULICS” 


by Enos J, Carlson and Carl R. Miller 
(Proc. Paper 953) 


ARTHUR I. MC CUTCHAN,! A.M. ASCE.—Under the heading “Measure- 
ment, Standard Samplers,” reference is made to the limitation of standard 
samplers in obtaining suspended material near the bed of a stream. In the 
authors’ reference 2, it is pointed out that the US D-49 sampler has a limita- 
tion at the other end of the range, in that it is basically unsuitable for samp- 
ling streams more than about 20 feet deep with two-way sampling, or 40 feet 
with one way sampling. Where rivers may in flood reach a greater depth 
than this, the need for sampling at these high stages is accentuated by the 
fact that, as shown very clearly in the authors’ figure 4, sediment concentra- 
tion tends to increase with discharge, so that the sediment load carried dur- 
ing flood is an even higher percentage of the annual load than is the flood dis- 
charge compared with annual stream flow, i.e. if it is important in stream 
flow studies to have accurate measurement of flood discharge, it is generally 
even more important in sediment studies to measure sediment load during 
fioods. 

This depth limitation of the US D-49 sampler is the main difficulty that 
has arisen in the small amount of work that has been done in Queensland. 

On the Dawson River where sampling was attempted in February, 1956, with a 
D-49 sampler modified to take Imperial pint milk bottles, the river rose toa 
near-record height of 72 feet, and in 12 days the stream discharge was about 
1,800,000 acre feet compared with the mean annual discharge of about 400,000 
acre feet. 

In an endeavour to overcome this difficulty without going to the complica- 
tion of the US P-46 sampler, an attempt was made to use a flexible plastic 
Jag inside a perforated rigid case which could be used in the modified D-49 
sampler. This approach is suggested in the authors’ reference 2, but no de- 
ails are given. The sampler, which had been obtained on loan from the Snowy 
Mountains Hydro-electric Authority, had to be returned before tests could be 
-ompleted, and the work was carried out only in a small laboratory flume, 
vithout any field testing. However, from a comparison of volume of sample 
ybtained with the flexible container and with the milk bottle under the same 
onditions, it appeared that satisfactory results could be obtained in stream 
relocities not less than about 4 feet per second, Presumably the decline in 
‘olume obtained in the plastic container at lower velocities is due to the in- 

; easing proportion of the velocity head of the water that is used in opening 
he bag. Best results were obtained with the air outlet plugged. The material 
ised for the bag was a poly-ethy! plastic 2 inch thick. 

_ The main advantage of this type of container is that it enables point samp- 
ing to be carried out by attaching to the nozzle of a D-49 sampler some sim- 
le type of mechanism as illustrated in the authors’ figure 2A, although the 


Executive Engr., Irrig. and Water Supply Comm., Queensland, Australia. 
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operation of the flap would be more difficult than on the hand sampler. Sud- 
den inrush effects are avoided without the complication of special provision 
for balancing pressures at all depths as is required with the rigid container. 

Reference to the results of research along these lines should be of con- 
siderable interest. 


SAM SHULITS,! M. ASCE. —The authors present an excellent summary of 
what is known and what should be known in sediment hydraulics. The need for 
short-range quantitative prediction methods, say for 50 to 100 years, is em- 
phasized. As the authors demonstrate, this is the realm in which the hydrau- 
lic engineer or the sedimentation specialist has achieved such remarkable 
progress- since the middle thirties when E, W. Lane and the writer in the 
Bureau of Reclamation ventured to predict the morphologic effects of Hoover 
Dam (then the Boulder Dam) on the Colorado River. | 

The twelve “phases of sediment hydraulics wherein research is needed,” — 
summarized by the authors on pages 953-18 to 19, point to other relevant 
research urgently needed, and described below. 


Fact-Finding Surveys 

The topic, “bed material transport computations,” is such a corner stone 
to the whole field of sediment hydraulics, that it is surprising that no au- 
thoritative and definitive comparative study exists to delimit the utility of the 
numerous formulas for the ped-load or the unmeasured sediment transport. 
Though two new formulas or computation procedures at least, have been 
added in the last few years, no one has come forth to defend, condemn or re- 
late the 14 formulas available at the last count of the writer. Experienced 
and capable sediment hydraulicians agree neither on the best formulas nor on 
the limitations of those used. A courageous fact-finding survey, with sup- 
ported conclusions, would certainly fill a practical research need and give 
the hydraulic engineer reliable knowledge on which formula to use when, 
where, and under what conditions. 


Fluvial Morphology 
The four items, “stable canals and channels, design for stability, aggrada- 
tion and degradation, and deposition in reservoirs,” extremely important in 
themselves and the subject of many current investigations, are also facets or 
components of the grand theme, fluvial morphology. Bank revetments will be 
local and temporary makeshifts, either in creeks or large rivers, if the mi- 
gration of the so-called “unstable” bend or bank is not considered. Perhaps — 
the alleged “instability” is partly due to the tortuous lateral movement of the 
whole stream. The aggradation and degradation which follow cutoffs or chan- 
nel rectification are part of the form adjustment of the profile caused by the 
drastic corrective measure, The morphology of these horizontal and vertical 
shifts of the stream is in need of considerable research, so that the overall 
plan and profile configurations can be predicted and correlated with the local 
corrective measures. ~ 
That such broad relationships can be found, is attested by the bold and q 
imaginative pioneering studies of the U. S. Geological Survey, in particular i 
those of Leopold and Maddock,(28) Wolman,(A)* and Leopold and Miller. ey 
— 4 
1. Civ. Eng. Dept., Pennsylvania State Univ., University Park, Pa. qi 
* To avoid confusion in the reference numbering scheme, the writer’s addi- : 
tions are marked A, B, etc. 7 
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Other important work on the quantitative analysis of watershed geomor- 
phology is that of Strahler.(C,D) Wittman,(E) after studying 100 years of re- 
cords for 250 miles of the Rhine River between Basle and Bingen, found a 
linear relation, M = KQ, wherein M is the migration of the crossings or bars 
in meters, K a coefficient, and Q the total discharge in cubic meters for the 
period under consideration. Inglis(F) has assembled empirical equations 
which relate river width, and the period and amplitude of river meanders. 
Strikingly systematic relationships for river profiles have been found by 
Yatsu(G) and the writer. (H,I) 


Bulk versus Particle Research 

The foregoing paragraph emphasizes “bulk” research, the sort in which 
broad and synoptic, yet quantitative laws are sought. This might be contrast- 
ed with “particle” research, (so-named for want of a better term), in which 
the goal or hope is to follow the behavior of the sediment particles as they 
are entrained into the suspended and bed load, dunes, anti-dunes, bars, banks 
and meanders. There exists a tendency to regard the latter class of research 
as “fundamental” or “basic,” and even more scientific. Actually there is 
need for both, as the two lanes of inquiry will certainly merge eventually. 
With our present limited fund of design criteria, the writer believes that 
“bulk” research should receive more attention. 


Deposition in Reservoirs 
Fortunately, the distribution of sediments in reservoirs is under study at 
oresent by a Task Force of the Hydraulics Division of the Society. 


Terminology 

A study of nomenclature in this field may soon be necessary. The authors 
introduce “sediment hydraulics,” Leopold and Maddock use “hydraulic geome- 
Ty,” while the writer likes “fluvial morphology.” The term “bed-load” is 
awkward, as it is used to designate both the transported material and the 
juantity. 

The U. S. Bureau of Reclamation is to be complimented on its work and 
ogress in “sediment hydraulics,” so ably described by the authors who 
hhemselves have been associated with these accomplishments. 
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Discussion of 
“TRANSITION PROFILES IN NON-UNIFORM CHANNELS” 


by Francis F, Escoffier 
(Proc. Paper 1006) 


ACHILLE LAZARD,!—The author’s paper is extremely interesting be- 
cause it permits the systematic extension of such interesting notions as those 
of “transition profiles” and “transition discharges.” 

The writer has, personally, no objection to the substitution of the word 
“transition” for that of “characteristic” which he had used in Ref, 1 of the 
paper. 

In creating the notion of “characteristic depth” Mouret(2) confined him- 
self to specifying that this depth is a characteristic of the channel bed and is 
ndependent of the discharge in contrast to the notions of normal and critical 
lepths which do depend on the discharge. 

In creating the notions of a characteristic discharge, which is the dis- 
-tharge for which the critical slope is equal to the slope of the channel, the 
yriter wished to bring out the fact that this discharge separates two types of 
low in the channel. The expression “transition discharge” is certainly to be 
referred, 

At the writer’s request the author has indicated that there can exist 0, 1, 
ry 2 transition profiles, depending on the slope of the channel. 

As the water-surface profile is normally horizontal at its point of inter- 
ection with the transition profile it appears to the writer indispensable to 
five a name or an index to each section of the water-surface profile that can 
ye defined. Mouret numbered the sections of water-surface profiles that he 
iad studied with the arabic numerals 1 to 6. (These correspond to the types 
3, Mg, My, S3, S2, and Sj, in the order named, as used among American 
mgineers.) The writer designated by the capital letters A, B, C, and C' the 
ew sections which he had introduced. These are shown in the accompanying 
igure. The sections C and C' are practically horizontal because of their 
ong inflection and appear to the writer to represent something real. 

_ In combining certain of Mouret’s sections with those of the writer, one ob- 
4ins water-surface profiles very similar to those that are described in the 
lassical manuals. However, certain of Mouret’s sections do not seem to the 
riter to represent real water surfaces. 

_ The author has in turn introduced new sections of water-surface profiles, 
ome of which probably do not represent real water surfaces. It would be 
seful to designate the others in some convenient way. 

The writer wishes finally to indicate that it is important in numerical cal- 
lations not to omit from the formulas the classical coefficient @ which ac- 
Bots for the unequal distribution of velocity about the mean velocity repre- 
ented by V. 

The coefficient a, about which a great deal is not known, should in princi- 
vary when the cross sections and the channel slopes vary. The evaluation 
this term can present difficulties. 


Ingenieur en Chef des Ponts et ChausSees, Ingénieur en Chef, French 
Railways, Paris, France. 
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In the case treated by the writer, where the channel is uniform and the 
channel slope is constant, a constant value could be ascribed to the coeffi- 
cienta. 

The introduction of this coefficient into the formulas modifies profoundly 
the numerical values obtained for the transition depths and discharges as 
these are determined by the intersections of lines that are nearly parallel. 
This is shown in the following table which contains the values calculated by 
the writer for the trapezoidal tailrace canal for the SOULOM powerhouse 


(Pyrenees), which has a slope of 0.004, for the extreme values for & of 1 anc 
1.2. 


Transition Transition 
Depth Discharge 
os _(in cm) _ {in m3/sec) 
1.0 90 7 
1.2 45 2. 


In closing, the writer will express the wish that the laboratories interest 
themselves in these problems and seek, by means of model experiments witl 
canals of appropriate cross sections, to verify or disprove the theory of 
water-surface profiles for gradually varied flow to which the author has 
made such an important contribution. 
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